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Abstract 

These are lectures presented at the Les Houches Summer School "Topology and Geometry in 
Physics", July 1998. They provide a simple introduction to non perturbative methods of field theory 
in 1 + 1 dimensions, and their application to the study of strongly correlated condensed matter 
problems - in particular quantum impurity problems. The level is moderately advanced, and takes 
the student all the way to the most recent progress in the field: many exercises and additional 
references are provided. 

In the first part, I give a sketchy introduction to conformal Held theory. I then explain how 
boundary conformal invariance can be used to classify and study low energy, strong coupling fixed 
points in quantum impurity problems. In the second part, I discuss quantum integrability from 
the point of view of perturbed conformal field theory, with a special emphasis on the recent ideas 
of massless scattering. I then explain how these ideas allow the computation of (experimentally 
measurable) transport properties in cross-over regimes. The case of edge states tunneling in the 
fractional quantum Hall effect is used throughout the lectures as an example of application. 
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Introduction 



Quantum impurity problems have been for many years, and increasingly so recently, a favorite subject 
of investigations, for theorists and experimentalists alike. There are many good reasons for that. 

First, these problems often represent the simplest setting in which some qualitatively fascinating 
physical properties can be observed and probed. For instance, the Kondo model (for a comprehensive 
review^] , consult |Q) provides a clear cut example of asymptotic freedom. The basic experimental fact is 
that normal metals with dilute impurities exhibit an unusual minimum in the temperature dependence 
of the electrical resistivity. The ultimate explanation is that the interaction of the electrons with the 
impurity spins produces an increase of the resistivity as T is lowered, counteracting the usual decrease 
in resistivity arising from interactions with lattice phonons. Low temperature means low energy, or 
large distance: we thus have a problem where interactions increase at large distance, a characteristic of 
asymptotic freedom Q As another example I would like to mention the recent experiments |2|, |j| about 
point contact tunneling in the fractional quantum Hall effect (for a review on this active topic, see Q]) 
at filling fraction v = i. Measurments of the shot noise reveal a behaviour (/ 2 ) oc |/g in the limit 
of weak backscattering, where few quasiparticles tunnel, and do so independently. If one compares this 
formula with the standard Schottky formula for Fermi liquids, one sees that this noise has to be due to 
the tunneling of fractional charges e* = |: although the existence of these (Laughlin quasiparticles) had 
been conjectured for a long time, the noise is the first direct evidence of their existence [J As a final 
example, let me recall that the basic archetype of dissipative quantum mechanics (for a review, look at 
||, Q), the two state problem coupled to a bath of oscillators with Ohmic dissipation, is described by 
another quantum impurity problem: the anisotropic Kondo model. Crucial fundamental issues are at 
stake here, as well as a large array of applications in chemistry and biology. 

The second reason of our fascination for quantum impurity problems is that they are, to a large 
extent, manageable by analytic methods. This has led to incredibly fruitful progress in the past. For 
instance, the renormalization group was, to a large extent, borne out of the efforts of Kondo, Anderson 
and Wilson to understand the low temperature behaviour of the Kondo model (see |Q). Also, the works 
of Andrei , Wiegmann || and others showed that the Bethe ansatz could be used to analyze situations 
of experimental relevance: this spurred a new interest in quantum integrable models, an area which, 
together with its various off-springs like quantum groups, knot theory and others, has become one of 
most lively in mathematical physics. 

Finally, it is fair to say that quantum impurity problems are not only of fundamental interest: they 
are at the center of the most challenging problems of today's condensed matter, like Kondo lattices, heavy 
fcrmions, and, maybe, high T c superconductors. 

In these lectures, I will concentrate on what is usually considered the most important about quantum 
impurity problems: their properties as strongly interacting systems. There is no doubt that strongly 
correlated electrons are of the highest interest. On the practical side, besides high T c superconductivity, 
let me mention the remarkable recent developments in manufacturing and understanding small systems 
like quantum wires, carbon nanotubes and the like, where, because of the reduced dimensionality, Fermi 
liquid theory is not applicable, and the interactions have to be taken into account non perturbatively. On 
a more fundamental level, strongly interacting systems exhibit rather counter intuitive properties, the 
most spectacular being probably spin charge separation. It is a challenge for the theorist to understand 
these properties, and quantum impurity problems no doubt provide the best theoretical and experimental 
laboratory to do so. 

Maybe it is time now to define what I mean by a quantum impurity problem. The general class of 
systems I have in mind have the following features: (i) There are extended gapless (critical) quantum 

1 These being only lecture notes, I have tried to refer to papers that were pedagogically inclined, if at all possible, rather 
than to original works. 

2 The analogy with QCD can be made more complete, including the logarithmic dependences encountered in both 
problems, and the "dimensional transmutation". 

3 The conductance itself is not a measure of the charge of the carriers. 
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mechanical degrees of freedom, which live in an infinite spatial volume, the "bulk" (ii) These interact 
with an impurity, localized at one point in position space. This impurity may carry quantum mechanical 
degrees of freedom. 

To have an example in mind, consider the Kondo problem: (i) The extended degrees of freedom are 
those of the bulk metal. The presence of a Fermi surface means that the metal sits at a RG fixed point 
(see e.g. ||). Physically this is easily understood since the system of electrons has ( "particle- hole" ) 
excitations of arbitrarily low energy about its Fermi-sea ground state, providing the critical degrees of 
freedom in the bulk of the metal (ii) The impurity spin, located at one point in space (say the origin), is 
a dynamical quantum mechanical degree of freedom (the dynamical process is the spin- flip). 

In the Kondo model we also see another feature of quantum impurity problems: they are generically 
one-dimensional. The problem of an impurity in a Luttinger liquid to be discussed below is inherently 
one-dimensional, but the Kondo model needs to be reduced to one dimension. Since the impurity spin 
sits only at one point in space, it is only the s-wave wavefunctions of the metal electrons that can interact 
with the spin. Second-quantizing this s-wave theory, we get a (non-interacting) quantum field theory 
of one-dimensional Fermions, defined on half- infinite (radial) position space (a half- infinite line), which 
interacts with the quantum spin at the end of the line. Considering a path-integral representation of the 
ID theory of Fermions, we have a (1 + l)-dimensional Lagrangian field theory, one dimension from the 
(half-infinite) radial space coordinate, and another dimension from the (say, euclidean) time coordinate. 
All interactions take place at one point in space, the end of the line, where the impurity is located. In 
the (1 + 1) dimensional space-time picture, the impurity sits at the "boundary" of space-time, which 
can be viewed as the upper-complex plane, the "boundary" being the real axis (these points of view are 
schematically illustrated in figure 1 and 2). This will be the general picture that we use in the study of 
all quantum impurity problems. 



Figure 1: We will think of quantum impurity problems in various ways. The figure on the left represents 
a f + 1 quantum point of view where the bulk right and left degrees of freedom are confined to a half 
line, with the impurity at the origin. Altrnativcly, because the theory is massless in the bulk, one can 
unfold this picture to get only right degrees of freedom on the full line, as indicated on the right. 




Figure 2: Alternatively again, one can go to imaginary time and obtain a 2D statistical mechanics point 
of view, with a theory defined in a half plane, that often we will chose to be the upper complex plane. 
In this figure, the arrows are supposed to "represent" the bulk degrees of freedom; later, we will see that 
they can be associated with integrable quasiparticles. 

In the Kondo problem we see one way how a quantum impurity problem can be realized experimentally: 
a bulk system (here the 3D metal) contains a finite but small concentration x of quantum impurities. 
In the limit of very dilute impurities (x < 1) the impurities do not interact with each other (to lowest 
order in x), and the single- impurity theory may be used to describe the physics of the bulk material in 
the presence of dilute impurities. Actually, experiments performed at very low concentrations are known 
to be in good agreement with the single-impurity theory for the ordinary (one-channel) Kondo model. 

A quite different realization of quantum impurity models occurs in the context of point contacts. 
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These are basically electronic devices: two leads (capable of transporting electrical current) are attached 
to a single quantum impurity. Each one of the leads is connected to a battery, so that electrical current 
is driven through the quantum impurity. One can then measure experimentally the electrical current J, 
flowing through the quantum impurity, as a function of the applied driving voltage V (from the battery) 
fl(i|| . The I(V) curve, the differential (non-linear) conductance Gdiff = dI(V)/dV, or the tempera- 
ture dependence of the linear response conductance G\ r = limy^ are examples of experimental 
probes characteristic of the quantum impurity. Notice that since the I{V) curve as well as the differential 
conductance are non-equilibrium properties, the point-contact realizations of quantum impurity prob- 
lems are theoretically more challenging than other realizations, in that more than equilibrium statistical 
mechanics is involved to achieve a theoretical understanding of these quantities. 

Ideally, the most interesting point contact situation would involve one-dimensional leads, where elec- 
trons are described by the Luttinger model, the simplest non-fermi- liquid metals jll| . It consists of left- 
and right-moving gapless excitations at the two fermi points in an interacting 1-dimensional electron gas. 
In the past, this model had been difficult to realize experimentally however. This is simply because in 
a one-dimensional conductor (such as a quasi-one-dimensional quantum wire so thin that the transverse 
modes are frozen out at low temperature), random impurities occur in the fabrication. These impurities 
lead to localization due to backscattering processes between the excitations at the two fermi points. In 
other words, the random impurities generate a mass gap for the fermions. 

Fortunately, there is another possiblity: the edge excitations at the boundary of samples prepared in 
a fractional quantum Hall state should be extremely clean realizations of the Luttinger non-fermi liquids, 
as was observed by Wen p2| . In contrast to quantum wires, these are stable systems because for 1/v an 
odd integer, the excitations only move in one direction on a given edge. Since the right and left edges 
are far apart from each other, backscattering processes due to random impurities in the bulk cannot 
localize those extended edge states. Moreover, the Luttinger interaction parameter is universally related 
to the filling fraction v of the quantum Hall state in the bulk sample by a topological argument based on 
the underlying Chern-Simons theory, and does therefore not renormalize. The edge states should thus 
provide an extremely clean experimental realization of the Luttinger model. 

We now describe an experimental set-up in more details ||l3| , Q (see figure 3). 



Figure 3: A schematic experimental set up to study point contact tunneling in the fractional quantum 
Hall effect (the magnetic field points towards the reader). Details are provided in the text. 

A fractional quantum Hall state with filling fraction v = 1/3 is prepared in the bulk of a quantum Hall 
bar which is long in the ^-direction and short in the y-direction. This means that the bulk quantum Hall 
state is prepared in a Hall insulator state (longitudinal conductivity a xx — 0), and that the (bulk) Hall 
resistivity is on the v — 1/3 plateau where a xy = (l/3)e 2 //i. This is achieved by adjusting the applied 
magnetic field, perpendicular to the plane of the bar. Since the plateau is broad, the applied magnetic 
field can be varied over a significant range without affecting the filling of v = 1/3. Next, a gate voltage 
V g is applied perpendicular to the long side of the bar, i.e. in the y direction at x = 0. This has the 
effect of bringing the right and left moving edges close to each other near x = 0, forming a point contact. 
Away from the contact there is no backscattering (i.e. no tunneling of charge carriers) because the edges 
are widely separated, but now charge carriers can hop from one edge to the other at the point contact. 

The left-moving (upper) edge of the Hall bar can now be connected to battery on the right such that 
the charge carriers are injected into the left-moving lead of the Hall bar with an equilibrium thermal 
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distribution at chemical potential Similarly, the right-moving carriers (propagating in the lower 
edge) are injected from the left, with a thermal distribution at chemical potential ^,r. The difference 
of chemical potentials of the injected charge carriers is the driving voltage V = hr — If V > 0, 
there are more carriers injected from the left than from the right, and a "source-drain" current flows 
from the left to the right, along the a;-direction of the Hall bar. In the absence of the point contact, the 
driving voltage places the right and left edges at different potentials (in the y-direction, perpendicular 
to the current flow), implying that the ratio of source-drain current to the driving voltage V is the Hall 
conductance G — ve 2 jh (both in linear response and at finite driving voltage V). When the point-contact 
interaction is included, at finite driving voltage, more of the right moving carriers injected from the left 
are backscattered than those injected from the right, resulting in a loss of charge carriers from the source- 
drain current. In this case we write the total source-drain current as I(V) = Iq(V) + Ib(V), where Ib(V) 
is the (negative) backscattering current, quantifying the loss of current due to backscattering at the point 
contact. It is this backscattering current that I ultimately want to show how to compute. 

Let me write up some formulas as a preamble. I will not have time in these lectures to discuss 
bosonization or edge states in the fractional quantum Hall effect: I will thus simply claim that, in its 
bosonized form, the problem is described by the hamiltonian 

h = ^J cfe [n 2 + (a^) 2 ] +Acos-^ L -0 i? )(o), (i) 



H \ \ 

where = v. Here, the free boson part describes the massless edge states |12J, and the cosine term 
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describes the effect of the gate voltage, with X oc V g . In general of course, the backscattering term 
induced by this gate voltage should be represented by a complicated interaction; but we keep only the 
most relevant term (the only one for v = i), which is all that matters in the scaling regime (see below) 
we will be interested in. 

The interaction is a relevant term, that is, in a renormalization group transformation, one has, b being 
the rescaling factor Q 

^ = (1-^)A + 0(A 3 ) (2) 

This means that at large gate voltage, or, equivalently, at small temperature (since then, elementary 
excitations have low energies, so the barrier appears big to them), the point contact will essentially split 
the system in half, and no current will flow through 0. The questions the theorist wants to answer are: 
how do we study the vicinity of the weak-backscattering limit? How do we find out more precisely what 
the strong back-scattering limit looks like? How about its vicinity? Finally, can we be more ambitious 
and compute say the current at any temperature, voltage and gate voltage? 

For this latter question, let me stress that we are interested in the universal, or scaling, regime, which 
is the only case where things will not depend in an complicated way on the microscopic details of the 
gate and other experimental parameters. In practice, what the experimentalist will do is first sweep 
through values of the gate voltage, the conductance signal showing a number of resonance peaks, which 
sharpen as the temperature is lowered. These resonance peaks occur for particular values V g — V* of 
the gate voltage, due to tunneling through localized states in the vicinity of the point contact. Ideally, 
on resonance, the source-drain conductance is equal to the Hall conductance without point contact, 
i.e. G resonance = ve 2 /h. This value is independent of temperature, on resonance. Now, measuring for 
instance the linear response conductance as a function of the gate voltage near the resonance, i.e. as a 
function of SV g = V g — V* , at a number of different temperatures T, one gets resonance curves, one for 
each temperature. These peak at 5V g = 0. Finally, these conductance curves should collapse, in the limit 
of very small T and 6V g , onto a single universal curve when plotted as a function of 5V g /T 1 ~ v ' . This is 
what the field theorist wants to compute. 

4 In particle physics language, 4r = —/3(A), so our relevant operator corresponds to a negative beta-function, ie an 
asymptotically free theory. 

5 This feature is actually remarkable. What it means is that, for one dimensional electrons with short distance repulsive 
interactions, an arbitrarily small impurity leads to no transmittance at T = p[lo[ : compare with the effects of barriers on 
non-interacting electrons you studied in first year quantum mechanics. 
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To proceed, it is useful to formulate the problem as a boundary problem. For this, a few manipulations 
are needed . We decompose = 4>l + 4>r and set | : 

<f e (x + 1) = -j= [<j> L {x, t) + M-x, t)} 

<p°(x + t) = -j=[4> L (x,t)-<t> R (-x,t)] (3) 

Observe that these two fields are left movers. We now fold the system by setting, for x < : 

(t>l = V e (x + t), 4> e R = ^{-x + t), 

4>i=<P°(x + t), cf>° R = -<p°(-x + t), (4) 

and introduce new fields 4> e '° = <f> e £° + (j> e jf , both defined on the half infinite line x < 0. The odd field 
(f> e simply obeys Dirichlct boundary conditions at the origin </>°(0) = 0, and decouples from the problem. 
The field <jf , which we call rather <!> in the following, has a non trivial dynamics 

H^H* = \J° dx [n 2 + a$)2] +A cos^. (5) 

The aspects we have to understand are, by increasing order of complexity: the fixed points, their 
vicinity, and what is in between. This is the order I will follow in these lectures. 



6 More details are given in Part I. A common objection to the following manipulations is that they are good only for 
free fields, but not when there is a boundary interaction. This, in fact, depends on what one means by "fields" - the safest 
attitude is to imagine one does perturbation theory in A. Then, all the quantities are evaluated within the free theory, on 
which one can legitimately do all the foldings, left right decompositions, etc. 
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Part I 

Conformal field theory and fixed points 

The first difficulty one encounters in that field is how to describe the low energy fixed points. This may 
sound rather simple in the tunneling problem, but in other cases, for instance in a tunneling problem 
for electrons with spin, the matter is quite involved. The reason for this is, that fixed points are not 
necessarily described in terms of nice linear boundary conditions for the bulk degrees of freedom. It does 
seem to be true however, that even if the quantum impurity has internal degrees of freedom, interaction 
and renormalization effects do turn the dynamical quantum impurity into a boundary condition on the 
extended bulk degrees of freedom, at large distances, low energy or low temperature. At low temperatures 
the system may be in the strong coupling regime (for instance, this is where Kondo's result diverges). 
The boundary condition is thus a way to think about the strongly interacting system. Nozieres' physical 
picture of screening |TJ| illustrates how this works for the simplest case, the one-channel Kondo model: 
the antiferromagnetic interaction of the impurity spin with the spin of the conduction electrons, which 
has renormalized to large values at low temperatures, causes complete screening of the impurity spin. 
A modified boundary condition on the electrons that are not involved in the screening, is left. This 
mechanism, however, appears to be much more general, and seems to apply to all quantum impurity 
problems. 



The boundary conditions generated in this process may be highly non-trivial (see e.g. [J17, 18 1). How- 
ever, since the bulk is massless (critical), the induced boundary condition is scale invariant asymptotically 
at large distances and low temperatures. Actually, it is, in most cases, conformally invariant. 

Quantum impurity problems are thus intimately related with scale-invariant boundary conditions: 
these are RG fixed points, and, like in bulk 1 + 1 quantum field theories, (recall that the bulk is always 
critical in the type of systems that we are considering here), conformal symmetry is the best way to 
describe them. 

Now, conformal invariance is a long story. All I can do is provide, in the next sections, what I believe 
is the minimal set of ideas necessary to understand what is going on, and tackle without fear the literature 
on the subject. In several instances, I will have to discard entire discussions of key issues, substituting 
them with some intuitive comments, and only providing the final result. Additional bits and pieces are 
then provided in the text in small characters, together with specific references, to help the reader bridge 
the gaps. Good reviews on this subject are the Les Houches Lectures of 1988 [[ljj], the article by J. Cardy 
pp| , the lectures by J. Polchinski pl[ , and the textbook ]2^] . The relevant chapters in |?3| can also be 
quite useful. 

In the following, I will intimately mix path integral and hamiltonian points of view. The two are of 
course equivalent, but each has its own advantages. 



1 Some notions of conformal field theory 
1.1 The free boson via path integrals 



We consider the free bosonic theory, with action 

S = — I dx\dx2 



(di$) 2 + (9 2 $) 2 



(G) 



To start, let us discuss briefly the issue of correlators and regularization. To keep in the spirit of condensed 
matter, we initially define the Gaussian model on a discrete periodic square lattice of constant a by setting 



Slo 



*(nO]' 



(7) 



O'fe) 



G 



where the sum is taken over all pairs of nearest neighbours. Introduce the lattice Green function 

1 ^ pik.r 

G latt (r) = jzY, 4 _2cos^-2cos^ ( 8 ) 

ni,7l2 L 

where the sum is restricted to the first Brillouin zone |rij| < and the prime means the zero mode is 
excluded [|. One has then (where the points r, r' belong to the lattice) 

(HrMr')) latt = G latt (r-r>) (9) 

while Gi a tt satisfies the discrete equation (where A; att is the discrete Laplacian) 

AlattGlatt = Slatt{r) (10) 

The important points here are the behaviours G; a tt(0) ~ — ^ In -j, L >> a, and Gi a tt{r) ps — In for 
a << r << L. 

We recall now Wick's theorem, according to which the average of any quantity can be obtained as a 
sum of all pairwise contractions. It follows that 



(e^e^^) iatt = exp [_ + ^? ($ 2 M) ia J 

x exp [-frh ((HrMr')) latt - M> Zatt )] (H) 



To define a continuum limit for this model, we look at distances large compared to the lattice spacing 
but small compared to L, where the right hand side of the previous expression simplifies into 

e ^) e ^(0\ w (ax<* + A)V<« /Jr-^*"* (1JJ) 



Latt \LJ 

The well known observation follows that the correlator vanishes unless charge neutrality is satisfied, that 
is (3\ + @2 = 0. We then have, where r, r' are now arbitrary points in the continuum, 

e **M e -tf*(r'>)„ [j-^rf^ (is) 

In the following we will sometimes, but not always, set a = 1. 
1.2 Normal ordering and OPE 

We now introduce complex coordinates z = Xi+ix2, z — xi—ix 2 . We have 8 = h{8\ — 182), 8 = ^(81+182), 
d 2 z = 2dx\dx2, and we define the delta function by J d 2 zS 2 (z,z) — 1, so S 2 (z,z) — ^5 2 (x\,X2)- The 
action reads now 

S= [ d 2 zd<$>8$, (14) 



and the laplacian 

A = 288 + 288 (15) 

Of crucial importance is the result [] 



8[ i ) =2ttS 2 (z,z) (16) 



7 The zero mode divergence simply occurs because the action is invariant under the symmetry <E> — » $ + est. 

8 A physical way to prove this is to observe that, in practice, — has to be regulated with a short distance cut-off which 

does introduce a z dependence, as for instance in H ( zz ~ a ) [ j{ t ne Heavyside function. 
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from which it follows thaiQ 

dd\n\z\ 2 = ddlnz = 2ir6 2 (z,z) (17) 
It is customary to write the basic correlator as0 

<$(z,z)$(z',z')) =-^\n\z-zf (18) 

Note that, in this expression, we have completely discarded the L dependence that occurs in the lattice 
system. A reason for doing so is that $ is not a "good" field anyway, and that we will usually consider 
rather derivatives of for which this ambiguity does not matter. The L dependence is however crucial 
for exponentials of the field $. When we discard it, we have to remember that, at the end of the day, 
only correlators which have vanishing charge are non zero. 
Now, still using Wick's theorem it follows that 

dB ($(z, z') ..)=- l -{5 2 {z-z'rz- z') ...), (19) 

where the dots stand for any other insertions in the path integrals, that involve no field either at z or z' . 
A relation that holds in this sense is simply rewritten 

<9<9$(z, z)$(z', z ) = -\s 2 {z -z',z- z'). (20) 

This is the first example of equations we will write quite often between "operators" in the theory - the 
word operator here occurs naturally when one splits open the path integral to obtain a hamiltonian 
description, see later. 

Recall that the equations of motion for the field $ read, on the other hand 

dd^(z,z)=0 (21) 

It follows that the product $(z, z)$(z', z') obeys the equation of motion except at coincident points. 

In the sequel, we shall constantly use the concept of normal ordering. We define the normal ordering 
of the product of two bosonic fields byp] 

: <J>(z, z)$(z', z') := <f>(z, z)$(z', z') + — In \z - z'\ 2 (22) 

47T 

This definition is such that the normal ordered product of fields now does satisfy the equation of motion 
even at coincident points 

dd:<f>(z,z)<!>(z',z'):=0 (23) 

As a result of this, the normal product is (locally) the sum of an analytic and antianalytic function, and 
can be expanded in powers of z. Thus, for instance 

$(z, z)$(0, 0) = — - In \z\ 2 + : $ 2 (0, 0) : +z : <9$$(0, 0) : +z : 9$$(0, 0) : + . . . (24) 

4-7T 

This is the first example of an operator product expansion (OPE). Like equation (j2fj), its precise meaning 
is that it holds once inserted inside a correlation function. OPEs in conformal field theories are not 
asymptotic, but rather convergent expansions; their radius of convergence is given by the distance to 
the nearest other operator in the correlation functions of interest. The right hand side of (|24|) involves 

9 Note that the two derivatives cannot be interchanged on singular functions, that is why Alnr = 2ttS 2 (xi, X2), and not 
twice as much. 

10 Of course the notation is somewhat redundant, since the value of z determines x\ and x^\ but in what follows, we will 
reserve the notation f(z) for analytic functions. 

11 If one wishes to keep the a and L factors, the normal ordering formula reads : <l?(z, 2)$(z', z') := <E>(z, z)<&(z' , z') + 



i In (\z ~ z'\/o) 2 . One has then (: $ 2 :\ = lni/a. 
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products of fields at coincident points, which turn out to be well defined in this theory. Notice that 
: <I> 2 (0,0) : in (|24| ) could be defined equally well by point splitting, as will be discussed later. 

For a product of more than two bosonic operators, the definition of normal order can be extended 
iteratively 

: $(z,z)$(z 1 ,z 1 ).. .<$>(z n ,z n ) := $(z,z) : $(z u z x ) . . . <$>(z n , z n ) : 

+ — (in \z — zi\ 2 : $(z 2 , z 2 ) ■ ■ ■ &(z n , z n ) : + permutations) (25) 

such that the classical equations of motion are still satisfied (a quick way to understand normal order, 
as is clear from the previous formula, is that quantities inside double dots are not contracted with one 
another when one computes correlators). 

Of crucial importance is the normal ordered exponential : e l/3 * :. It is a good exercise to recover the 
OPE 

. e i/3i$( z> 2) . . e »/3 2 *(0,0) ._ . e i/3i*( Z ,z)+i/3 2 $(0 : 0) . ^6) 

and 

d$(z, z) : e i/3 *(°'°) := ^- : e^°> > : + : d$(z, z)e 1 ^ ^ : (27) 

The quantities inside the normal ordering symbols can now be expanded in powers of z, z like an ordinary 
function. 

Exercise: Show that the "quantum Pythagoras" theorem holds: 



cos2Vtt < !' :) 2 + (: sin2^/iT^ 



finite 



Notice that in (|27|), we could have treated 9$ as an analytic function: a perfectly legitimate thing 
to do when one computes correlators. Hidden here is, in fact, the very convenient decomposition of the 
field itself into the sum of an analytic and antianalytic component (one has to be very careful however 
when one writes $(^, z) — 4>(z) + 4>(z); first, because such a decomposition does not hold for the general 
fields summed over in the path integral, and second, because the field <& does not obey the equations of 
motion at coincident points). Pushing that line of thought a bit further however, one has 

(4>{z)4>{z')) = -^Hz-z') 

{^zW(z>)) = --L\n(z-z>) (28) 

this up to phases due to the branch cuts. We will also use in the following the dual field 

$ = 0-0 (29) 

The exponentials are scalar operators, ie they are invariant under rotations. More general operators 
also have a spin, so their two point function reads 

<0(*>*)°(°.o)> = (30) 

The simplest example is provided by O = 9<J> which has h = l,h = 0. In general, single valuedness of 
physical correlators requires h — h to be an integer. The numbers h and h are called usually right and 
left conformal dimensions. The dimension of O is d = h + h, while its spin is s = h — h. 
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1.3 The stress energy tensor 

The stress energy tensor is defined in the classical theory as follows. Consider a coordinate transformation 
— ► Xfi + £/j (that is, changing the arguments of the fields in the action from x^ — > x^ + e M ) . The 
variation of the action reads, to lowest order ^|, 

5S = —-^ J d l _ L e v T^ y dx 1 dx 2 (31) 

Elementary calculation shows that 

V - -lird^dv® + nS^dpQdp® (32) 

The stress energy tensor in the quantum theory is defined through Ward identities: the end result 
is the same formula as ([32]), but where products of fields are normal ordered. It enjoys some very 
important properties: the symmetry T\% = T%\ as a result of rotational invariance, and the tracelessness 
Tn + T22 = as a result of scale invariance. In addition, the stress energy tensor is always conserved, ie 
the operator equation d^T^ = holds. This can be checked explicitely for the free boson using ( |32| ) and 
the classical equations of motion, which we recall are satisfied by normal ordered products^]. In general, 
one introduces complex components 

T zz = \(T u -T 22 -2iT 12 ) 
Tez = \ (Tn - T 22 + 2zTi 2 ) 
T Z - Z = T- ZZ = i(T u +T 22 ). (33) 

which, in the free boson case, read simply 

T zz = -2tt : {84>f : (34) 
T n = -2tt : (84>) 2 : (35) 
T z - Z = 0. (36) 

(37) 

Using the equations of motion, T zz is analytic (again, in the special sense that, when inserted in correlation 
functions, the dependence is analytic away from the arguments of the other operators) and will simply 
be denoted T(z) in what follows; similarly, T zz is antianalytic. It is easy to chek that these properties 
extend to other models: a remarkable consequence of locality (so there is a stress energy tensor to begin 
with) and masslessness is thus the existence of a conserved current, a field of dimensions (2,0). Currents 
provide a powerful tool to classify the fields of the theory, as we will see shortly. 

Plugging back our results into (|3l|), one checks that the variation of the action SS vanishes exactly 
for a conformal transformation z — ► z + e(z). This is the celebrated conformal invariance, which we will 
discuss in more details below. 

The short distance expansion of T with itself reads 



T(z)T(0) = + 4r(0) + -9T(0) + analytic (38) 

2z 4 z z z 



The coefficient of the l/2z 4 term is uniquely determined once the normalization of T has been chosen. 
For other massless relativistic field theories, this coefficient takes the value c/2 where c is a number known 
as the central charge. For a free boson we see that c = 1. For n independent free bosons, c = n. For a free 
Majorana fermion c=\. The relative normalization of the two other factors is fixed by the requirement 
that T(z)T(0) = T(0)T(z). 

12 The factor of 2n is peculiar to the conformal field theory literature. 

13 More generally, that T is classically conserved follows simply from the fact that the action is stationary as a consequence 
of the classical equations of motion, so 8S must vanish for arbitrary e. 
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Exercise: show this. 

The fact that T is analytic except when its argument coincides with the one of some other field inside a correlator has 
an interesting consequence for the trace of the stress energy tensor. Indeed, d(T(z)T(0)) is now a derivative of the delta 



function! Using the conservation equation d^T^, which reads in complex components, &T - 
that the two point function of the trace of the stress energy tensor is non zero: 



([T 11 +T 2 2}(x 1 ,x 2 )[T 11 



-T 22 ](0,0)> = — AS 2 (x u x 2 ) 



jc9(Til + T 2 2) = 0, if follows 



(39) 



This is a simple example of an anomaly, a . q uantity which is zero classically, but non zero quantum mechanically. It is a bit 
dangerous to give too much meaning to however, since the trace is not really an independent object - it is much safer 
to remember again that classical equations of motion hold except at coincident points. 

Remarkably, the stress energy tensor was introduced in a statistical mechanics long ago by Kadanoff 
and Ceva |2^] . These authors were interested in the way Ising correlators change under shear and scaling 
transformations. They recognized that, in a critical theory, rescaling in the x\ and X2 directions was 
equivalent to changing the horizontal and vertical couplings, and thus that the effect of shear and scaling 
could be taken into account by introducing an operator "conjugated" to these changes in the correlators, 
just like, say a change in temperature could be taken into account by introducing the total energy in the 
correlators, ft is thus possible to physically identify T, and to wonder how its continuum limit behaves, 
and how the various algebraic properties we are going to derive emerge. 

1.4 Conformal in(co) variance 

To fix ideas, let us now consider the free bosonic theory on a cylinder of circumference y and length L. 
Introducing the complex coordinate w = x + iy such that the imaginary axis is parallel to the cylinder's 
length (sec figure 4), the two point function of the field <£> in that geometry is easily found to be 



4-7T 



sin7rT(u; 



w') 



(40) 



x2 



w = x + i y 



Figure 4: Some of the geometries used in the text. 
From this, it follows that 



cylinder \ | sill 7rT(w — Ul') \ 



/3 2 /2ir 



Let us focus on the holomorphic part 



. e WM . . e -i04>(w') ■ 



ttT 



cylinder \ Sm nT(w — W 1 ) 



/3 2 /4tt 



This can be shown to be, equivalently, 

cylinder 



e i/3(j>(w) . . e -i04>(w') 



dz\ h ( dz'^ 



plane 



(41) 



(42) 



(43) 
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where we used the mapping z = e — 2i7rTl ". The latter formula expresses the covariance of the two point 
function under the conformal transformation. Another example of such covariance is provided by the 
derivative of the field 



(d w <f>(w)d w ,<f>{w')) = (|0 f|0 (d z <j>(z)d z ,<f>(z')) , (44) 

where we suppressed mention of the geometry, which is implicit in the variables used. 

Relations like (43 4^) are well expected, since the gaussian action is, in fact, conformal invariant; 



this follows, as discussed above, from the properties of the stress energy tensor, and thus is expected to 
generalize to other local massless field theories. More directly, this invariance is easily established for the 
free boson, since upon changing the argument of the field from z — > w in the action, S is invariant, the 
Jacobian cancelling the term coming from the partial derivatives. 

Of course, one has to be quite careful in using the conformal invariance of the action, since the 
correlators are not invariant - ie, one has for instance, (&(w, w)Q(w', w')) ^ (&(z, z)Q(z', z')), while the 
naive change of variables in the action would suggest the propagators to map straightforwardly, and thus 
the equality to hold. The reason for this discrepancy comes from the cut-off, which is also modified in 
a conformal transformation. We, on the other hand, wish to use the same regularization whatever the 
geometry, ie, for instance, use a square lattice of constant a to regularize both the problem in the plane 
and on the cylinder; hence, there is an "anomaly" . 

Fields obeying the general covariance relation (and a similar one for the antiholomorphic part) 

(0(w)0(w')) = (2L\ f£L\ (0(z)0(z')) (45) 

are called primary fields. The field <f> itself is not primary, though, in a way, it satisfies the equivalent of 
the previous relation with h = 0, since 

(Hw)Hw')) = {<j>{z)<t>W)) + ^ In (^^) (46) 

Fields which are not primary exhibit in general more complicated covariance relations. An example is 
provided by the second derivative of 0, which we leave to the reader to work out. A more interesting 
example is furnished by the stress energy tensor. Though we have defined it so far by normal ordering, 
it is clear that an equally good definition is obtained by point splitting, ie rj 



T(z) = -2?r lim 

d-,0 



We have thus, from the change of variables 



d<j>{z + d/2)d<f>(z - d/2) 



And 2 



(47) 



T(z) = -2tt lim 



w'(z + d/2)w'(z - d/2)d w <j){w(z + d/2))d<j>(w(z - d/2)) 



And 2 



To define the stress energy tensor on the cylinder, we use the same definition (47), but with z replaced 
by w 0. Therefore 

T(z)^[w'(z)] 2 T(w) + ^{w 1 z}, (48) 



where the added term is 



1 1 w'(z + d/2)w'(z - d/2) 1 

12 {WlZ > ~ d™ 2 [w{z + d/2) - w(z - d/2)} 2 ~ 2^2 



14 Here, it does not matter how and by what amount the two points are split of course, provided they both tend to z at 
the end 

15 That is, normal ordering is always defined by subtracting the short distance, geometry independent, divergences. 
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This is known under the name of Schwartzian derivative, and reads 

. . 2w"'v/-3w' 4 

{w,z} = — 49 

It enjoys nice properties under the composition of successive conformal transformations, that we leave 
to the reader to investigate. An important property following from ( ]48| ) is that T acquires a finite 
expectation value on the cylinder, while it did not have one in the plane 

{T) cylinder = ~^ T f (50) 

Exercise: show this by using the propagators on the plane and the cylinder, together with appropriate 
definitions of normal ordering. 

The OPE of the stress energy tensor with any field of the theory has the general form 

T(z)O(0, 0) = . . . + Ao(o, 0) + -dO(0, 0) + . . . (51) 

Z* Z 

The two terms explicitely written follow from the fact that O has dimension h and the use of the Ward 
identity (|57|) . It can be shown that 1/z 2 is the highest singularity if O is primary. 
Exercise: check this for the free boson by considering various examples. 

1.5 Some remarks on Ward identities in QFT. 

Suppose in general that there is a transformation of the field (xi,x%) = &(xi,X2) + S$>(xi,X2) that leaves the product 
of the path integral measure and the Boltzmann weight e~ s invariant. Examples of such transformations are provided for 
instance by translations or rotations in ordinary isotropic homogeneous physical systems. Consider then a transformation 
<E>'(x'l , X2) = &(xi,X2) + p(x\ , X2)&Q(x\ , 2:2). For general p, this is not a symmetry of the problem anymore. On the other 
hand, we can always change variables in the functional integral and reevaluate any correlator in terms of the new field <&'. 
This means that we have the identity 

0=1 [d#']e" s ' - / [d<S>]e- s (52) 



where in S' one maybe had to add up terms coming from the change of variables in the path integral. On the other hand, 
one can expand the right hand side of this equation to first order in the change of fields assumed small. Sincefor p a 
constant the product of the measure and the weight would be invariant, this means that the right hand side of ( fed ) must 
depend on the gradient of p only, ie one has 

rhs= f[d$]e~ s J j^d^ pdx 1 dx 2 (53) 

The quantity is called a Noether current. Since it comes from local manipulations, it must be a local quantity. Now, 
that (B3J) vanishes is something that must hold for any reasonably smooth function p. Let us choose p to be enual to unity 
inside a disk of radius Hi, and to vanish on and outside of R2, while it is arbitrary in between. Integrating ( Kg ) by parts, 
we get an integral of on the circle Hi, together with an integral on the annulus between Hi and H2 of d^j^. Since the 
functions p is quite arbitrary there, it follows that the current has to be conserved, that is 

dyh = (54) 

This result would still hold of course with fields inserted far from Hi and H2, so ( ]54[ ) truly holds as an operator equation, 
in the sense explained above. 

As an application, consider a translation x M — > x^ + e^, where is small and constant: we obtain a current which, in 
the classical case, coincides with ija = The foregoing procedure is a generalization to the quantum theory, and the 

conservation equation follows from Noether's theorem. 

Now consider some field O inside the circles, say right at the origin. Under the transformation <t> — > <E>', this field 
becomes O 1 = O + SO (the change is expanded to first order as before, but of course O' might as well depend on the 
derivatives of p in general). We now have, since all we are doing is changing variables in the path integral 

= / [d$']0' e~ s ' [d*]0 e~ s ... (55) 



Expand this to first order. This time, the integration by part gives 



iSO = ^- / (dx 2 ji - dxij 2 )O(0, 0) = Jr- I (dzj - dzj)O(0, 0) (56) 

Z7T I ZZ7T 
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As an application, consider a transformation 2^2 + a(z). where a is small. The Noethcr current associated with it 
is given by j = ia(2)T(2) and j = ia(z)T(z). For a transformation that is conformal inside a contour CL a nd diffcrcntiably 
connected to a (necessary non conformal) transformation vanishing at large distances, one finds from ( |5q ) the key result 

SO = — a(z)T(z)O(0)dz ~ — a(z)f(z)O{0)dz (57) 



2«7r J c 2iir Jc 

Exercise: derive from this ([u]) and the fact that there are no singularities stronger than I/2 2 for a primary operator. 

Notice finally that for the free boson, the expression of the stress energy tensor is almost the classical one, up to normal 
ordering, and it appears as if the integration measure essentially plays no role in the construction of the Ward identities. 
That one can forget about the behaviour of the measure in conformal transformations is justified a posteriori, by the fact 
that the quantum currents are indeed conserved. The measure would play a more subtle role for theories defined on curved 
two-dimensional manifolds. 



1.6 The Virasoro algebra: intuitive introduction 

As noticed before, the main consequence of conformal invariance is the existence of a conserved current, 
the stress energy tensor T. In general, one sets 

71 — — OO 

that is, plugging this expansion into the OPE provides a definition of what the field L n O actually is: for 
instance L O = hO, L_iO = dO, L 2 T = 1, etc. In general, one does not expect fields with negative 
dimensions to appear, or at least not fields with arbitrarily large negative dimensions (weird things can 
occur in non unitary theories adapted to disordered systems in particular though). This means that for 
every field, L n O must vanish for n positive large enough. Of particular interest are the primary fields, 
for which the highest singularity is 1/z 2 , ie they satisfy L n O = 0, n > 0. 

For the moment, we can contend ourselves with the intuitively reasonable notion that the L n are 
"operators" acting on the space of fields of the theory - here, exponentials multiplied by normal ordered 
polynomials in derivatives of the field $ - so the L n are not unlike differential operators acting on functions 
(in fact, the (L„i)™ are just that). 

It is then tempting to ask oneself what the algebra of these operators is, that is how do L n (L m O) 
compare with L m (L n O)? This is easily done by using contour integration together with the short distance 
expansions. The commutator [L n , L m )0 can be computed as follows. We have 

L n L m O = ^z n+1 T(z) J ^w m+1 T(w)O(0) (59) 

where the contours encircle the origin and C is inside C (see figure 5). 

Indeed, imagine writing the OPE of the integrand. First we expand T(w)O(0) to extract the field 
L m O, on which the action of L n is then obtained by the second integration. That the countour C is 
inside C is natural from the point of view of "radial quantization" which, as we will see later, gives a 
precise operatorial definition to the L„'s. It is also necessary if one wishes to use the OPE in the order we 
just said for convergence reasons. The product L m L n O is computed in the same fashion with this time 
with a contour C" inside C: this forces one to expand first the product T(z)0, resulting in the opposite 
order for the operators. Comparing the two, and forgetting the operator O itself, we see that 

[L n ,L m ]= [ [ ^-z n+1 w m+1 T(z)T(w) (60) 
Jc 2i7r Jc, 2l7r 

and a computation using the OPE of T with itself gives 

[L n , L m ] = (n - m)L n+m + j^i™ 3 ~ n)8 n+m (61) 

This is the celebrated Virasoro algebra, the L n being called Virasoro generators. It is an infinite dimen- 
sional Lie algebra. 
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Cz = C'-C 



C 




Figure 5: The contour manipulations that lead to the definition of a commutator in radial quantization. 

Exercise: compute the action of the Virasoro generators of say derivatives of the field cf>, and directly 
check the Virasoro commutation relations. 

A very important use of this algebra is to provide one with a natural structure to organize and 
recognize the fields in a theory. Of course, one does not quite need this powerful tool for the free boson, 
whose fields are easily built "by hand" , but for more complex theories, this is really very useful. Given 
a lattice model with microscopic variables, arbitrary combinations of neighbour variables can be built, 
whose scaling limit may or may not give rise to new scaling fields: which are truly new, which are nothing 
but "L„" 's (descendents) of others? What happens is that a theory has a certain number of primary 
fields, which is very often finite (eg, three for the Ising model), and all the other fields are just descendents 
of these ones. The whole set of fields is thus organized into products of representations of the left and 
right Virasoro algebras, for which the primary fields are heighest weight states. This can be expressed 
by the compact form 



The situation is quite similar to the case of angular momentum in ordinary quantum mechanics, where 
the space of say the possible electronic states of some atom can be organized in terms of representations 
of the angular momentum algebra. Of course, here we have an algebra with an infinite number of 
generators, instead of three for angular momentum in three space dimensions. Qualitatively, there is 
an infinite number of Virasoro generators because there are an infinite number of elementary conformal 
transformations, one for each power of z: z n . 

As in the theory of angular momentum, unitarity contrains the quantum numbers, that is the values 
of conformal weights for a given central charge. This in turn gives rise to strong constraints for multipoint 
correlation functions; this is beyond the scope of these lectures, but not by far. In particular, correlations 
at strong coupling fixed points in the Kondo model can be computed just by using that technique |Q. 

The reason why we focused on the algebra of the L n 's is because of the special role of the stress 
energy tensor in conformal transformations. Of course, we could define other operators and other algebras 
associated with any field that has an integer dimension (so the contour integrals can be closed in the 
complex plane. Generalizations also occur for fields with non integer, rational dimensions, and a cut 
plane, but this is more complicated). A natural candidate in condensed matter is provided by currents, 
for instance d(j>. Set therefore 




(62) 



h,h 




1 



1 



(63) 



From 



(dcj){z)d<f)(w)} 



1 1 



(64) 



An (z — w) 



2 ' 
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it follows that 



J rn-\-7i 



(65) 



Here we recognize the oscillator algebra standard in the quantization of the free boson, and maybe it is 
time to discuss more what the mode expansion has to do with hamiltonian quantization. 



1.7 Cylinders 

I shall mostly discuss what happens in the case of the cylinder. The key idea here is to remember that 
path integrals are scalar products of states. If we insert a field O at y = — oo on the cylinder, this 
corresponds to having prepared the system in a state (an "in" state) |0), to which is associated, in the 
field representation a wave function ^o^s 1 ]; result of a partial path integration 



*o[$ s i] = 



[d$n]0(-oo) e 



-Sn 



(66) 



3> s i fixed 



where the integral is taken over all configurations of the field in the bottom part of the cylinder f2 = 
(— oo, 0] x S 1 , the values at the boundary S 1 being held to $ s i 0. Similarly if we insert a field O' at oo, 
this corresponds to projecting the system on an "out" state \0'}, to which is associated a wave fucntion 
^O'i&s 1 )- The scalar product of these two states is then 



(OIO') 



d$si]*o'( $ si)*o(*si) ; 



(67) 



and this is essentially the correlation function of the two fields O, O' . Of course, by translation invariance 
on the cylinder, this does not depend on the particular place where we have cut open the path integral. 

To make things concrete, let us discuss an example we will use explicitly later, with O = O' = /, the 
identity operator - ie, nothing is actually inserted at ±oo. In this case, ([37]) is just the partition function 
Z of the problem. To find the wave functions, let us split open the path integral at y = 0, and let us 
Fourier decompose 

" (68) 



3v o) = $ " c 



where oj n — 2imT. Introduce then the solution of the Laplace equation A$o = subject to the constraint 
$o{x, y = 0) = <&(x). One finds easily 



(69) 



We can now split the field in the path integral into $ = $o + < &' where vanishes at y = 0. Because 
of this, together with the fact that <E>o solves Laplace equation, integration by parts shows that the path 
integral factorizes into the partition fucntions of two half cylinders with Dirichlct boundary conditions 
(we will get back to these later; the point here is that they are independent of ^s 1 ); an d an interesting 
term 



Z = Zl |[d$ sl ]exp Ll^|a; n ||$ n | 2 ^ 



Comparing with (67), it follows that, up to a phase 

ttT 2 



[&S 1 ] 0^ ex P 



where we Fourier transformed back the integrand in ( |70| 

16 1 am not being too careful here about what happens at — oo. 



$(x) - $(£') 

i \sin7rT(a; — x') 



(70) 



(71) 
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We notice here as a side remark that the issue of finding wave functions is more than formal: the 
computation above was carried out for instance by people interested in finding the wave function of the 
Thirring model in terms of the original fermions P5| . 

In this point of view, we have a Hilbert space made up of states in one to one correspondence with the 
various fields of the theory. For operators other than the identity, one has to be a little bit careful. While 
the in state is always obtained by inserting O at — oo, the out state is obtained actually by inserting 

(e-^) 2 * (e- 2i * T ™) 2 ~ h O at +00. 

The same analysis can be carried out in the plane in the framework of "radial quantization" , where 
time is In \z\. This is why the expansions we used earlier were called OPE. To be correct however, they do 
have a meaning as OPE's only when the operators are radially ordered, since recall that, to an Euclidian 
Green function computed with a path integral, there corresponds a time ordered Green function in the 
quantum field theory. Of course, other hamiltonian descriptions (for instance, the standard one where 
imaginary time runs say along X2) could be obtained by splitting open the path integrals differently. 

The remarkable thing now is, that the hamiltonian on the cylinder has a very simple spectrum. Indeed, 
first observe that for primary operators, using the formula (^3j), the two point function on the cylinder 
decreases at large distance along the cylinder as 

(O(w)O(w')) « exp [-271-TO + h){y - y')] cxp [-2inT(h - h)(x - x')] (72) 

For non primary operators, the same can be shown to hold, the more complicated terms in the covariance 
formula decreasing more quickly. 

On the other hand, suppose we want to describe the theory on the cylinder in a hamiltonian formalism 
with imaginary time along the y axis. As is well known, the rate of decay of correlation functions is given 
by the gaps of the hamiltonian, and the decay (|72| ) indicates that there is an eigenstate of the hamiltonian 
whose eigenenergy is 2irT(h + h) over the ground state. More generally, the computation of any physical 
property of the theory boils down to evaluating correlators, which all obey ([72]): therefore, the whole 
space of the quantum field theory must be organized in states associated with the various observables, 
such that their eigenenergy is 2irT(h + h) over the ground state! This is exactly what we expected from 
the hamiltonian formalism described before, with one additional piece of information: the spectrum of 
H. 

In addition, it is important to stress that the stress tensor acquires a non vanishing expectation value 
on the cylinder, due to the schwartzian derivative ( |50|) . As a result, the hamiltonian on the cylinder 
reads 

H = 2nT (l + L - (73) 

and the momentum 

P = 2-kT (Lq - Lo) (74) 

where Lq has eigenvalues h, Lq eigenvalues h. 

We can of course define the whole set of Virasoro generators on the cylinder by 

T(w) = (2ttT) 2 L n e- 2iwnTw - (75) 

Now we have a precise meaning to give the L n as operators, and their commutator can be computed, of 
course giving rise to the same Virasoro algebra derived more intuitively before. Notice the identity 



1 



l/T 



H=— (T + T)dx (76) 
27r Jo 

This is independent of result of analyticity. 

It should be clear that the whole conformal invariance analysis could be written within the hamiltonian 
formulation. For instance, the OPE of T with itself corresponds to the commutator 

-!- [T(x),T(x')\ = 6{x - x')T'(x) - 25' {x - x')T(x) + U"'(x - x'), (77) 

2l7T 6 
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and the OPE of T with a primary field 



— [T(x), 0{x')\ = 5{x - x')d x O - 6'(x - x')hO 

Z7T 

1.8 The free boson via hamiltonians 

We now discuss the hamiltonian formalism more specifically for the free boson. The Lagrangian is 



(78) 



L = l 
2 



l/T 



dx 



(d t $) 2 - (d x <s>) 2 



from which the momentum follows 
and the standard hamiltonian 



n = a$ 



h = 



i 

2 Jo 



l/T 



n 2 + (d x $y 



with the canonical equal time commutation relations 

[<f>(x,t),Il(x',t)} = iS{x~x') 



(79) 

(80) 
(81) 

(82) 



The field is periodic in the space direction, that is Q(x,t) = &(x + l/T, t). We chose to compactifiy the 
field on a circle of radius r, that is we identify $ = $ + 2irr. The mode expansion of the field reads then 
(see eg @ for many more details on this) 



$(x, t) = x + Tpt + 2irTrwx H = - ( a 



2iTrTn(x-t) _ - 2iTrTn(x+t) 
6 CX— n € 



(83) 



where x — T J $(x, t)dx — $o is the boson zero mode, while p = J n(x, t)dx = Ho is the total momentum. 
w in an integer (the winding number); p is quantized such that rp = k is also an integer. The commutation 
relations of the operators are [x,p] = i and 



[a n ,a m ] = [a n ,a m ] = nS n+m , [a n ,a m ] = 



(84) 



The operators a n ,a n are related with the usual creation and annihilation operators of the free boson 
harmonic oscillators by 

— 1 (85) 



and 



-i\/na n ,n > 0, a n = i\J —na_ n ,n < 



a n — ~i\/na- n ,n > 0, a„ — i\f —na\,n < 



(86) 



Note that if we go to euclidian space time, replacing t by — iy and then use the conformal coordinates 



2wT(y-ix) 



(z = e 



-2inTw 



„2irT(y+ix) 



), we obtain the expansion 



. , p wr\ . p wr 
$(z, z) = $n - i — H Inz - i — 



In z + — \= V - (a n z " + a n z n ) 



4?r "Tfn n 



(87) 



When the winding number is non zero, the field is not periodic around the origin; rather, a "vortex" is 
inserted there. Wh 
general, we will set 



inserted there. When w = 0, we can set ao = «o = 774^! one recovers the expansion ( |63| ) for id(j>. In 



P 



4tt 



P 



4tt 
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The hamiltonian (plf) reads , before regularization 



H = 2ttT 



, \2 P 2 1 I - \ 

n(wr) + h - > j (a- n a n + a-„a n ) 

47T A * — ' 



(89) 



A question that arises now is the relation between the normal ordering defined in the field theory 
and the normal ordering in the usual sense of ordering free bosonic operators in quadratic expressions: 
: a n a m :— tti n f(n,m) a sup(n,m)- The two might differ by a constant; in the present case actually, they 
coincide provided one uses zeta regularization. Indeed, by ordering H, we encounter a divergent term 
Yin, which we can regularize by (for results on the zeta function, see p7| ) 

f> = C(-l) = ~ (9°) 
i 

With this prescription, the hamiltonian with the vacuum energy divergence subtracted reads as it should 
(|73|), with the Virasoro generators 

m m 

together with 

1 x 1 v 

L o = ^a'o + ^ a -nCt n , L = -al + ^ "-n<5n (92) 

n—1 n—1 

The modes a n and a n are annihilation operators for n > and creation operators for n < 0. The whole 
space of fields is thus obtained starting from highest weight states \w,k) which are annihilated by the 
annihilation operators and are eigenstates of the zero modes, and applying creation operators to them. 
Schematically, one has 

H = Heis w ,k <8> Heis w< k (93) 

w,k 

Of course, the \w, k) states are primary, and thus highest weight of the Virasoro algebra. Accordingly, 
one could as well build the whole space of fields by acting on them with L' n s. This would be more 
complicated; for instance, the field dcj), which is simply the result of a_i |0,0), is not obtained from the 
action of L„i on that state at all. This means in general that more primary fields are necessary than the 
|w, k) in the Virasoro description. 

1.9 Modular invariance 

A convenient way of encoding the field content of the theory is to write the torus partition function, that 
is, the partition function when one imposes periodic boundary conditions in the imaginary time direction, 
too. One has, using (|73| ) 



Z = Tr exp 

Using the mode decomposition, one finds easily 



-2.71 (Lo + Lo- j£ 



(94) 



z = —^—-j2i hak ^ k ( 95 ) 

where q — e~ 27rTL — q (the notation allows consideration of more complicated parallelograms), 

oo 

V(q)=q^l[(l-q n ) (96) 
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and 

„ / k wr\ - ( k wr 

= 9.TT 1 h,..u = 9.TT 



h wk -2^— + -) ,h wk =2^—- T ) (97) 

An important property of this partition function is that it is modular invariant. What this means is, 
suppose one considers quantization of the free boson with time in the x instead of y direction. The radius 
being the same, this will lead to the same expression as (^5|) but with L and 1/T exchanged, that is 

Z= . * E(gO h " fc (gf" fc (98) 



where q' = e 2n / TL , The expressions (Bq) and (98) do turn out to be equal thanks to some elliptic 



functions identities (see next section) . They ought to be, of course, since they represent the same physical 
object from two different points of view. 

For more sophisticated theories, the partition function cannot be computed a priori, but it is possible 
to determine it by imposing that it does not depend on the description, ie is modular invariant. Sec p2| 
and references therein for more details. 



2 Conformal invariance analysis of quantum impurity fixed points 
2.1 Boundary conformal field theory 

An excellent reference for this part is the original work of J. Cardy (ps[]). Consider now a field theory 
defined only on the half plane x-i > (figure 6) - it might be for instance the continuum limit of a 2D 
statistical mechanics model which is at its critical point in the bulk, that is T = T c , the usual critical 
temperature of the system. Various situations could occur at the boundary depending on whether the 
coupling there is enhanced, or whether some quantum boundary degrees of freedom have been added. 



x 2 

n 




Figure 6: The geometry for boundary conformal field theory. 

Consider, to fix ideas, the simplest case where the statistical mechanics model would have the same 
couplings in the bulk and the boundary ( the so called "ordinary transition"). Intuitively, one expects the 
system to still be invariant under global rotations, dilations and translations that preserve the boundary, 
and that this invariance should be promoted to a local one, ie conformal invariance in the presence of the 
boundary. 

Physical fields are now characterized both by a bulk and a boundary anomalous dimension. If both 
fields are taken deep inside the system, they behave as in the bulk case. On the other hand, if they are 
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near the boundary, one has, for example, 

(0(xi,x 2 )0(x' 1 ,x' 2 )) w 1 - , \xi - x[\ » x 2 ,x' 2 (99) 

\xi — a; 1 |^ a " 

ie the large distance behaviour of the correlators parallel to the surface is determined by the boundary 
dimension. We quote also the formula 

(0(^X2)0(^,4)) « j^- (c 0S ef°- d (ioo) 

A condition of boundary conformal invariance is that T\ 2 = T 2 \ = when x 2 — 0, which means 
physically that there is no flux of energy through the boundary. As a result, the left and right components 
of the stress tensor are not independent anymore, but T = T for Imz — 0; this is expected, since the 
theory is invariant only under the transformations that preserve this boundary, that is satisfy w = w for 
Imz = 0. As a result however, one can define formally the stress tensor in the region Imz < by setting 

T{z) = f(z), Imz < (101) 

Instead of having a half plane with left and right movers, we can thus equivalently describe the problem 
with only right movers on the full plane. For instance, the two point correlation function in the half plane 
is related with the four point correlation function in the full plane. Also radial quantization corresponds 
to propagating outwards from the origin in the upper half plane, with hamiltonian (see figure 7) 

/ T{z)dz + cc (102) 
Jc 



x 2 
i I 




x 1 




Figure 7: Geometry of the contours for the boundary case. 

Using the continuation ( |l0l| ), this becomes a closed contour integral of T only: thus, the Hilbert space 
of the theory with boundary is described by a sum of representations of a single Virasoro algebra this 
time: 

H = J2 Vir h ( 103 ) 



The natural mapping in this problem is w = — — In z, which maps the half plane onto a strip of width 
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L with the same boundary conditions on both sides 0. The hamiltonian now reads 



H 



L 



(*-a) 



(104) 



Note that there are, roughly, two factors of two differing from the periodic hamiltonian: the prefactor 
has a 7r instead of 2n, and there is a single Virasoro generator in the bracket. The space onto which the 
periodic hamiltonian ([73|) acts is uniquely defined by the (bulk) theory one is dealing with, say the Ising 
model - as we discussed, this specification amounts to giving the various representations of Vir ® Vir 
defining the model. For (104), the space depends on the boundary conditions; it is specified by a set of 
representations of a single Virasoro algebra. By a careful study of the two point function in the plane, 
together with the conformal transformation (where the jacobians still involve the bulk dimension), one 
can show that the gaps of H are given by the corresponding surface dimensions. 

It is important to stress again that the same physical obervable will be associated with different 
representations of the Virasoro algebra in the bulk and boundary cases. For instance, the spin in the Ising 
model coresponds, in the bulk, to Vw*x/i6 <S> Vir 1 /i 6 , while with free boundary conditions, it corresponds 
to Vir 1/2 (with fixed boundary conditions, the spin is the same as the identiy operator). 



2.2 Partition functions and boundary states 

To classify boundary conditions, it is extremely useful to deal with partition functions a bit. We consider 
thus a cylinder with a periodic direction of length 1/T and a non periodic one of length L: on either 
side, boundary conditions of type a, b have been imposed. We can describe the situation in two ways 
(sec figure 8): cither imaginary time runs in the direction parallel to the boundary ("open channel"), in 
which case we can write the partition function as 



Tre 



-H ab /T 



(105) 



where H a b is the hamiltonian (104) with boundary conditions a and b, or imaginary time can run in the 
direction perpendicular to the boundary ( "closed channel" ) , in which case 



Z = {B a \e- LH \B b ) 

where \B a ) , \ Bb) are boundary states, and H is the periodic hamiltonian ([? 



(106) 



open 

4 



t 



Figure 8: The open and closed channel geometries when boundaries are present. 

Observe that the the bou ndary states are not normalized: they are entirely determined, including their 
norm, by the condition that ( 106 ) gives the right partition function. To make things more concrete, fixed 



17 Notice that here I have put L in the mapping, instead of T. In the periodic case, we could have used the mapping 
w = — 2j^r In z to produce a similar result. This is all equivalent, but I prefer the present choice, where 1/T is always the 
periodic direction in the problem. 

18 A strip with different boundary conditions on either side would correspond to a half plane with different boundary 
conditions x < and x > 0, with a "boundary conditions changing operator" inserted right at x = 0. 



22 



boundary conditions in the Ising model for instance are represented, in the microscopic Hilbert space, by 



the state LB, 



fixed 



Yli while for free boundary conditions one has \B) free = ]J i (|+) + |— )). 



Here the boundary states are states in the Hilbert space of the bulk theory, ie in Vir (g) Vir. Conformal 
invariance at the boundary requires 

(L n - L_ n ) \B) = (107) 
A solution to this equation is provided by so called Ishibashi states |3Q] 



\h) = ^2 \ h > n ) ® \ h ' n ) 



(108) 



where \h, n) denotes an orthonormal basis of the representation Virt, and |/i, n) the corresponding basis 
of Vim- 

In the case of the free boson, a boundary state will satisfy (107) if it satisfies a stronger constraint 

(a n ±a- m )\B) =0 (109) 

This in fact correspo nds t o Neumann and Dirichlet boundary conditions, for which T 12 oc c* 1 $92 ( f ) = 0. 
The negative sign in (10G) is solved by 



\B) oc exp 

Therefore, we can build boundary states by 



E 



E 



c k exp 



E 



a- n a- 



\0,k) 



\0,k) 



(110) 



(111) 



The question of interest is to determine the coefficients c&. A quick way to proceed is to recognize 
here a Diri chle t state: indeed, suppose we act with &(x,t = 0) on the boundary state. Because of the 
condition (109), the oscillator part just does not contribute; what does contribute is only the x part, 



which acts as x 



' dp 



Therefore, we have 



\b d {$ q ))=n d J2 

k=—oc 



e~ 4fc * o/r exp 



E 



Ct- n a- 



\0,k) 



(112) 



The last question, which is actually of key importance for what follows, is the determination of the 
overall factor Af: in other words, what is the overall normalization of boundary states? This is where the 
consideration of partition functions is useful. 

To answer this, we observe that, if we compute the partition function with height <&o on both sides, 
the identity representation should appear once and only once. On the other hand, the partition function 
is easily computed in the closed channel from the boundary states: one finds, for more general pair of 
values at the boundary 



-. oc 

Z = {B D {$ Q )\e- LH \B D {& Q ))=j^ D -— ]T q 



k 2 /8wr 2 p ki($ -$' )/r 



(113) 



k=— oo 



where q 



-A-kTL 



We now perform a modular transformation to reexpress this partition function in 



terms of the other parameter q = e 7r / LT . One has (the proof of this is a bit intricate. See eg |29| , 
chapter 3.) 

V (q) = ^=V(Q) (114) 



19 This topic goes back to the early days of open string theory. A nice recent paper on the subject is ]31| , where the 
following computations are carried out in many more details. 
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and, by using Poisson resummation formula for the infinite sum, 

exp (—iran 2 + bnj = —= cxp 



b 
2iir 



one finds 



i r (* -*0+ 27 ™'") 2 



(115) 



(116) 



This expression has a simple interpretation: one sums over all the sectors where the difference of heights 
between the two sides of the cylinder is $o — $o + 27rrn. For each such sector, the partition function is 
the product of a basic partition function corresponding to heights equal (without the 2irr identification) 
on both sides, times the exponential of a classical action. The latter is easily obtained: the classical field 



is $ 



$0 — <&o+27mr 



y, whose classical action is 



exp 



1 



2LT 



($o - $o + 2mir) 2 



Tr 



Vir h q 



-c/24 



Consider now (|116|). We know that the partition function must write as a sum of characters (that is, 



as follows from (103) and (104)) of the Virasoro algebra with integer coefficients; even 



though I will not spend time discussing what the characters at c = 1 are (q h /r) for generic h), it is 
easy to see that this implies that the prefactor in (116) has to be an integer. Since we do not expect 
the normalization of the boundary states to change discontinuously with <&o, this integer is actually a 
constant, whatever $ ,$g. We can in particular choose <£>o = 3>o, f° r which the identity representation 
Virh=o appears in the spectrum; of course it should appear only once, and therefore 



M D = 



1 



V27 



(117) 



The other condition corresponds to Neumann boundary conditions, or, equivalently, Dirichlet bound- 
ary conditions on the dual field $ = $o- One finds the boundary state 



e -2i-,rw<S> exp 



W — — QO 



E 



Ct-nOi- 



\w,0) 



The Neumann Neumann partition function reads then 

Znn 



and one has 



A/jv = x V 2 ' v 



(118) 



(119) 



(120) 



The Neumann Dirichlet partition function is actually independent of the values of &o, 'to, since then 
the field cannot wind in any direction. 

Exercise: show the following 

The consideration of boundary states is extremely powerful to find out and study boundary fixed 
points. A general strategy is, knowing the Virasoro algebra symmetry of the model at hand, to try 
to find out combinations of Ishibashi states that are acceptable boundary states. Solving this problem 
involves rather complicated constraints. For instance, if one has several possible candidates |-Bj), the 
partition function with boundary conditions ij can easily be evaluated in the closed channel; after modular 
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transformation to the open channel, it should expand as a sum of characters of the Virasoro algebra with 
integer coefficients. Another constraint is that the identity representation should appear at most once in 
all open channel partition functions. Clearly, this becomes a rather technical subject; more details can 
be found in the paper of J. Cardy J3^]. Questions like the completeness of boundary states (ie whether 
all the boundary fixed points of a given bulk problem are known) are still open in most cases. 

2.3 Boundary entropy 

Let us now suppose that we have a one dimensional quantum field theory defined on a segment of length 
L, with some boundary conditions at x — and x = L. As is well known, the partition function at 
temperature T of this theory will be given by the same expression as the partition function of the two 
dimensional systems considered previously; notice however that I have changed conventions calling now 
x (resp. y) what was y (resp. x) previously (see figure 9) 0. 

y = it ,. 



Figure 9: The geometry for defining boundary entropy. 

h 

In terms of the parameter q, this partition function is expressed as a sum of terms — with integer 
coefficients: the spectrum of H ab is discrete, and its ground state has integer degeneracy, nothing very 
exciting. In the limit L — > oo, the spectrum becomes gapless however, and one has to be more careful 
about the concept of degeneracy. If we take this limit, the free energy of the quantum field theory behaves 
as 

F=-Lf + f a + h (122) 

where / is a free energy per unit length, f a , fb are boundary contributions. These contributions will 
involve, as T — > 0, a boundary energy that is non universal, but also a boundary entropy. It is easy to 
see what this entropy will be by using a modular transformation. The same partition function expresses 
then as a sum of -4=r with some non integer coefficients that come from the Poisson resummation formula 

77(g) o 

(in general, from the modular S matrix). In the large L limit, q — > 0. From the fact that 

F = -TlnZ (123) 

we see that, as T — > 0, / = 0(T 2 ) (the ground state energy of H ab is set to zero in this approach; for the 
exact dependence of / on T 2 see section 6), while f a and f b are of the form 

f a = -Tlng a ,f b = -T\ng b (124) 

where ||: 

g a = (B a \0) , g b = (0\B b ) (125) 

A one dimensional massless quantum field theory defined on a line with boundary conditions (or boundary 
degrees of freedom as we will see next) therefore has a non trivial zero temperature boundary entropy, 
or ground state degeneracy. 

20 This is to match as much as possible with the literature; in any case, there is no perfect notation that would be 
convenient all the way through. 
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Exercise: Show that the precise meaning of this degeneracy is related with the behaviour of the density 
of states 




(126) 



where we parametrized the excitation energies of H a b by e n = n,L large (when computing the 
partition function and its logarithm, do not forget to integrate the fluctuations around the saddle point!). 



As we have seen in the previous subsection, some boundary conditions have a degeneracy g < 1, ie 
a negative boundary entropy. This is a bit shocking, but of cours e we should remember, first, that g is 
more a prefactor in an asymptotic formula for degeneracies (126) than a true ground state degeneracy 
(at L = oo, there is no gap), and second, that we are dealing with quantum field theories and that this is 
only a finite, properly regularized "entropy". The same remark applies, somehow, to g being non integer. 
However, it is perfectly possible to have non integer degeneracies for semi-classical systems involving 
kinks il. 
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Intermezzo 

Perturbation near the fixed points 



A scale-invariant boundary condition is a RG fixed point (recall that the bulk is always critical 
in the type of systems that we are considering here). As with any RG fixed point, there is a set of 
relevant/marginal/irrelevant boundary operators (and couplings) associated with each scale- invariant 
boundary condition. These operators have support only at the boundary, i.e. at one point in position 
space (at the position of the impurity). 

If no relevant boundary operators are allowed, then the scale invariant boundary condition represents 
a stable fixed point (the zero temperature fixed point, describing the Kondo model at strong coupling, 
is an example; so is the Dirichlet fixed point in the tunneling problem, to which we will get back soon). 
Irrelevant boundary operators give perturbatively calculable corrections to physical properties evaluated 
at the RG fixed point. Many important physical features of the Kondo model are actually due to the 
effect of the leading (dominant) irrelevant boundary operator ]l8[ . 

Adding a relevant boundary operator to the Hamiltonian describing a particular scale-invariant bound- 
ary condition, destroys that boundary condition, and causes crossover to a new, scale-invariant boundary 
condition at large distances and low temperatures (in the infrared). In other words, we have a (boundary) 
RG flow, describing the crossover from the initial scale-invariant boundary condition (in the ultraviolet, 
i.e. at short distances or high temperatures) to a new scale invariant boundary condition (in the infrared, 
i.e. at large distances and low temperature). 

Note that at every stage of this flow, the bulk remains always critical and unchanged; the only action 
is at the boundary. An interesting observation concerning general boundary RG flows was made in |33| : 
the zero-temperature boundary entropies (s = In g in the previous section) generally obey 

suv > sir, (decrease of boundary entropy) 

This may be viewed as a boundary analogue of the well known c-theorem of bulk conformal field theory 
|35| . (Note, however, that the universal numbers sjjv and sir do not seem to be obviously related to a 
dynamical quantity, in contrast with the central charge, which is related to the stress tensor of CFT). 

A well known example is the one-channel Kondo model. Initially, at weak coupling (at high tem- 
perature, in the ultraviolet), we have a quantum mechanical spin decoupled from the electron degrees 
of freedom of the metal. An isolated (s — 1/2) spin has a zero-temperature entropy of suv = hi 2. 
At strong coupling (at low temperature, in the infrared), this impurity spin is completely screened by 
the conduction electrons. This means that no dynamical degrees of freedom are left, and thus we have 
sir = 0. 

3 The boundary sine-Gordon model: general results 
3.1 The model and the flow 

We consider now the model we had decided to tackle in the introduction 

S=\J dxj dy [(3 X $) 2 + (d,$) 2 ] + A J di/oos|*(0,|/). (127) 

This model is called the boundary sine-Gordon model since it has a sine-Gordon type interaction, but at 
the boundary. In more general terms than those of the edge states tunneling, the physics of this model 
is rather clear. The limits A = and A = oo are fixed points, corresponding to conformal invariant 
boundary conditions, respectively of Neumann and Dirichlet types. Away from these limits, the model 
is not scale invariant because of the boundary interaction. In the vicinity of A = 0, the RG equation is 

^ = (l-.g)A + 0(A 3 ), (128) 
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where we have set g = v = j^. It is natural to expect that A flows all the way from to oo under 
renormalization. Equivalently, the boundary conditions look like Neumann at very high energy (UV) but 
like Dirichlct at low energy (IR) - the dimension of the physical coupling is [A] = L 9 ^ 1 , so the typical 
energy scale for the cross over between UV and IR behaviours is Tb oc A~. Equivalently also, the field 
$ feels Neumann boundary conditions close to the boundary, but feels Dirichlct boundary conditions 
instead far from it, with a cross over distance 1/Tg. 

Notice that the boundary entropies of the UV and IR fixed points are different. To compute them, 
we can use the results of the previous section after having identified the radius of the boson. In the shift 
$ — > $ + 2irr, the interaction cos must be unchanged, which requires 

r = \ (129) 



It follows that 



Notice the ratio 



- '■/ J / 02 \ 1/4 



2\ -1/2 

(131) 



9n_ = (P_ 
9d V 8?r 

For the case of a relevant perturbation we are considering here, this ratio is larger than one: the boundary 
entropy is greater in the UV than it is in the IR. This is in agreement with the intuitive idea that degrees 
of freedom disappear under the renormalization group, leading to a loss of information. There is a 
well known conjecture stating that for any allowed flow in a unitary system (that is, roughly, a system 
with real, local hamiltonian), g uv > gm- For the case of irrelevant perturbation, one finds g^ < go, so 
according to this the flow should not be possible, which is indeed the case: since the operator is irrelevant, 
it does not generate any flow, and one should observe N boundary conditions both at small and large 
distance. 



3.2 Perturbation near the UV fixed point 



The first question we will be interested in is the calculation of the boundary free energy at any temperature 
T and coupling A. This can be represented by a Coulomb gas expansion as follows. First, by using a 
conformal mapping, one finds the two point function of the free boson with Neumann boundary conditions 
on the half cylinder 

sin irT(y - y') 



(132) 



<<%)<%')) = -9ln 

Exercise: derive this, by first computing the two point function on the half plane. 

We can then evaluate the ratio of partition functions with and without boundary interaction as follows 

Z{\) — i r^/T 

rT = 1 + 

n=0 



Z{\ = 0) 



1 + E(2^! A2n / d yi ...dy 2n {cos ^( yi )... cos ^(y 2n )y (133) 



Of course, only electrically neutral configurations with n positive and n negative charges contribute. 
After some rescaling, one finds 



Z(X = 0) 



71=1 



where the dimensionless coupling is 



2T 



(134) 



(135) 
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and the integrals are 




Figure 10: Charges of a two dimensional Coulomb gas that move on a circle d = 2 sin " 2 2 Ul . 

The integrands have small distance behaviours l/u 2g . It follows that there is no short distance 
divergence, and the integrals are all finite for g < | (there are never large distance divergences here 
since we have a temperature). When g > |, the integrals have divergences. In the sequel, I will always 
regularize integrals dimensionally, not by introducing a cut-off. To explain what this means, consider the 
case n = 1, which can be done by elementary computations 

T r(i - 2g) 

/2 - ■ (137) 

This can then be continued beyond g — \ simply by using the known continuation of T to negative 
arguments. How to do this in the case of arbitrary n is a bit more tricky. A way to do it relies on 
the remarkable fact that the integrals lm can be expressed in an almost closed form by appealing to 



techniques of Jack polynomials 1 36 , B7l . I will only give the result here 



1 " / y \rm + g(n — i + 1 )] \ 2 

hn = W ? H ( v rK + ff (n- l ) + i]J (138) 

where the sum is over all sets (Young tableaux) m = {m\, . . . ,m n ) with integers m.j obeying mi > 
TO2 . . . > m n > 0. This expression can be used to compute the l2n numerically to high values of n, or, 
more fundamentally, to perform the analytic continuation in g. I will not discuss this further, and get 
back now to the physics of this model. 

For g < i at least, the perturbative expansion is well defined, giving a series in A with positive 
coefficients. This series will presumably have a finite radius of convergence - although one does not 
expect the appearance of a singularity on the positive real axis (this would correspond to the existence of 
a phase transition on the one dimensional boundary). Beyond this radius, some other technique has to 
be used to understand quantitatively what happens. It is possible to argue what the leading behaviour 
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of the partition function at large A should be. Indeed, Z actually depends only on the ratio of the two 
energy scales Tb and T, so large A is like small temperature. But small temperature corresponds, going 
back to an euclidean description, to a cylinder of large diameter. In this limit, the partition function 
per unit length of the boundary should have a well defined, "thermodynamic" limit, so Z should go as 
Z oc exp (^r-). This means our perturbative series has to go as expicstX 1 / 1 ^ 9 ). 



3.3 Perturbation near the IR fixed point 



A natural idea to find out what happens beyond the radius of convergence is to think of the problem 
from a "dual" point of view, ie around the A = oo infra red fixed point. The first question one may ask 
is along which irrelevant operator this fixed point is approached. There are several, equally interesting 
ways to answer this question. 

The first one starts by considering the model where the bulk degrees of freedom have been integrated 
out, leading to the action (at zero temperature, which makes the formulas a bit simpler) 



2tt 



dydy' 



y-y' 



+ A 



J dy cos ^(y) + j J dy{d y $>f 



(139) 



where we have added an irrelevant mass term to make some integrals finite. 

We can find kinks interpolating between adjacent vacua and satisfying the equations of motion [|38| 



m<9,,$ = sin — $ 

y 2 2 



A simple solution of this equation is indeed 



„ , >. 27T 8 _i 

fins{y) = -j + -j tan 



exp - 




(140) 



(141) 



The energy of this kink is infinite, but can be made finite by subtracting a constant term from the 
action, replacing the cos ^$ by cos ^$ — 1. If we then consider a configuration of the field $ made of a 
superposition of far apart instantons and anti-instantons, 

* = I>/*».(v-Wi) (142) 

the kinetic term of the action can be conveniently evaluated by Fourier transform 



Shin 



/ l*HI 2 l-l^ 



At large distances, one finds 



Ski 



16tt 



j<k 



(143) 



(144) 



This is in exact correspondence with the Coulomb gas expansion discussed previously, but with the 
exchange ^ — > It follows that the IR action reads, at leading order 



S 



dx / dy 



(145) 



where we recall that $ is the dual of the boson <E>. It follows that the IR fixed point is approached along 
an operator of dimension h = ^ = |. One also checks that A^ oc A~ » . 

It is important to stress now that, while the flow away from the UV fixed point is fully specified by 
a single perturbing term, the situation is very different for the approach towards the IR fixed point. Of 
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course, one is f ree i f one wishes to perturb the D boundary conditions by a single irrelevant operator as 
represented in ( 145 ), though of course one has to be especially careful in defining the theory because of 
the strong short distance divergences in the integrals. The point is, that there is only one particular way 
of approaching the IR fixed point that corresponds to the trajectory originating at our UV fixed point. 
This means that the large A behaviour of the series we are interested in would be computable from the 
knowledge of an action of the form 



dx / dy 



(dy$ 



_\ 2 



47T ~ 

dy cos— $(0,y) 



k 



A" 



(146) 



where Ok belong to a very large class of operators allowed by symmetry: there are for instance all the 
cos^$, (d<&j , and many others. Since all these operators come with appropriately scaled powers of 

the coupling constant, they all give contributions to physical properties that depend on ou r sin gle scaling 
variable A, and no operator can be discarded (let me stress that an expansion such as (146) does not 
make much sense until one specifies the regularization procedure employed) . 

Here of course the reader should ask: but why didn't we add that collection of operators near the 
UV fixed point as well? The point is that we had control of what we wanted to do near the UV fixed 
point, and only a maniac would want to use such an irrealistically finely tuned combination of operators 
to perturb a fixed point. However we have no control about the way the IR fixed point is approached: 
this is entirely determined by the dynamics of the quantum field theory, and it turns out to be quite 
complicated. It is important in particular to realize that, starting from (146) and trying to go against 
the renormalization group flow, there is, most probably, only one choice of IR perturbation that would 
get back to our UV fixed point. 

It seems very hard to push the instanton expansion beyond the first non trivial order to try to get at 



(146), or get the expansion of physical quantities for large A: this has for a long time made IR perturbation 
theory impossible to carry out beyond the first trivial order. 

Remarkably however, (146) can be entirely determined using ideas of integrability. In a scheme where 
everything is dimensionally regularized, the only vertex operator that is present near the IR fixed point 
is cos ^$: none of the other harmonics actually appear! There are also very strong constraints on the 
other operators. 

In any case, the non-perturbative region of large A is very hard to access quantitatively using pertur- 
bation of the IR fixed point. Fortunately, the problem can be tackled by using ideas of integrability, a 
topic to which we will turn soon. 



3.4 An alternative to the instanton expansion: the conformal invariance anal- 
ysis 

Clearly, the instanton expansion, if physically appealing, is a bit laborious, especially when one considers 
how little information it finally provides. The conformal invariance analysis gives an alternative way, 
usually more reliable, to know which operators are present near the IR fixed point. Indeed, this informa- 



tion is encoded in the partition function Zod (116) for <I>o = ^' : the modes in the sum correspond to 



operators with dimension h = ™ , the cosines identified previously, while the other terms obtained by 
expanding the eta function correspond to powers of derivatives of the field. 
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Part II 

Integrability and the complete flow 

The constraint of integrability has been used with much success to study crossover scaling in bulk 2D 
theories exactly. However, it is often objected that integrable models are not so relevant for experimentally 
observable physics for at least two reasons: 

(a) : In order to achieve integrability, extensive fine-tuning of parameters is often required. Therefore, it 
is often believed that exact predictions made by studying an exactly integrable model might often not 
be generic and therefore difficult to observe experimentally. For example integrable spin chains with spin 
greater than 1/2 are gapless, while the generic even-spin spin chain and those observed have gaps in the 
spectrum. 

(b) : A very important set of experimentally accessible observables are transport properties. Amongst 
those is the conductance which is usually computed from (equilibrium) Green's functions using the 
Kubo-formula. It is usually very difficult or impossible to compute exact Green's functions, even when 
the system is (Bcthe-Ansatz) integrable. Therefore, before the progress made in the last couple of years, 
integrability was largely restricted to the computation of thermodynamic quantities, excluding transport 
properties at finite temperature. 

In fact, the situation is quite different for integrable quantum impurity problems. For these problems, 
exact transport properties can be computed (even out of equilibrium) and integrability can answer directly 
experimentally important questions. 
In particular: 

(i) : In order to achieve integrability in quantum impurity problems, one often needs to adjust very few 
parameters - sometimes none! For instance, both the Kondo effect and the point contacts in fractional 
quantum Hall Effect devices, provide an experimental realization of an integrable system without any 
fine-tuning. Similarly, while integrable higher spin quantum spin chains are non generic, the higher spin 
Kondo problem, is integrable. 

(ii) : Exact transport properties (at non-zero temperature) can be computed from integrability |f39f . The 
linear response conductance for the quantum Hall point contact, for instance, agrees quantitatively with 
recent experiments by Milliken/Webb/Umbach [[l3): this is the quantity we will discuss in what follows. 
Many other properties can also be computed, and I will discuss them briefly at the end. 

The method to obtain those exact results is a bit unconventional: it relies crucially on a judicious 
choice of basis of the Hilbert space of the system. We use a basis that is natural from the point of view of 
integrability. It is simply the basis in which all the infinite number of conservation laws (that exist since 
we have an integrable system) are diagonal. This basis turns out to have a "Fock-space" like structure, 
i.e. it is spanned by "quasiparticles" . It is in this basis that the quantum impurity interaction becomes 
tractable. In order to compute transport properties, we use a kinetic equation for those quasiparticles of 
the Bethe-Ansatz. This is non-trivial, since we are really describing a fully interacting system, where a 
single-particle (Fermi-liquid) concept such as a kinetic equation seems out of place, at first sight. The 
single particle kinetic equation would fail to produce exact results in interacting systems due to the 
existence of particle production processes in the single-particle basis. However, the quasiparticle basis 
dictated to us by integrability is precisely characterized by the absence of quasiparticle production (and 
"factorized scattering"). This particular and special feature of an integrable theory, allows us to use a 
kinetic equation to compute transport properties exactly from integrability. 

A historical note is necessary here. The point of view I will use is different from the original works on 
integrable quantum impurity problems. In the latter works, the authors started from a "bare" theory, and 
proved by hand that there were simple eigenstates obtained by making a two body scattering ansatz, the 
Bethe ansatz (there won't be much about the Bethe ansatz per se in these lectures) . They then proceeded 
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to build the "physical" theory by filling up the ground state, and studying excitations above it. Nowadays, 
it has become customary to start directly with the physical theory and prove its integrability using a 
very different and powerful tool that I will introduce below, perturbed conformal theory. The spectrum of 
excitations and the S matrix are then deduced (some would say guessed) by using symmetry arguments; in 
particular by analyzing non local currents (for this aspect, see eg f40|| ). The approach gives rise somehow 
naturally to the computation of transport properties, in particular by making physical sense of massless 
scattering, and that has definitely been a progress. Another key advantage is that, for a given "physical" 
theory, there are many possible "bare" choices, that is many different possible regularizations. Usually 
only one of them is integrable, and not always the obvious one; for instance, S. Zamolodchikov showed in 
his pioneering work that the Ising model at T c with a magnetic field is an integrable quantum field theory, 
but it is well known that its standard regularized square lattice version is not integrable (in the case we 
are interested in, the boundary sine-Gordon model, I am actually not aware of any simple integrable 
bare hamiltonian). All this is not to diminish the beauty and astonishing insight of the pioneering works 
about quantum impurity problems ^ ^| , nor the huge body of work on the Bethe ansatz and Yang Baxter 
equation that it is impossible to even start to acknowledge here. 



4 Search for integrability: classical analysis 

As emphasized in the first sections, one of the main uses of conformal invariance is to provide a convenient 
basis to the Hilbert space of observables in terms of representations of the infinite dimensional symmetry. 
For the free boson, the basis furnished by irreducible representations of the Virasoro algebra is just one 
of many choices: the basis furnished instead by representations of the Heisenberg algebra is also possible, 
and sometimes more convenient. 

When one wishes to study the problem with a boundary interaction, the question arises, of which 
basis will be the most convenient to work. It turns out it is still a third choice, provided by a "massless 
scattering theory" . To understand what this means, it is good to first consider the classical case. 

In the classical limit, one can scale the parameter j3 off the action. Going to real time, one obtains a 
classical scalar field $(x,t) satisfying the Klein-Gordon equation in the bulk x S [— oo,0): 

<9 4 2 $ - d 2 $ = (147) 

together with the boundary conditions (where A is also rescaled): 

(148) 

Now, the point is that the "most natural" basis of solutions for the bulk problem, that is plane waves, 
behaves badly with respect to the boundary interaction: if a plane wave is sent towards the boundary, 
what bounces back is a complicated superposition. Is it possible to find a better basis made of wave 
packets that will bounce nicely? 

To find such a basis, we make a detour through a more complicated problem which everybody knows 
is integrable, the massive sine-Gordon model. That is, we want to think of our Klein-Gordon problem as 
the A — > limit of the sine-Gordon equation: 

<9 4 2 $ - <9 2 $ = -Asin($) (149) 

It turns out that this massive model, in the presence of the boundary interaction, is still integrable. We 
will discuss this point in more details below, and start instead by considering how things look like in the 
massless limit. 

There are two types of finite-energy solutions of the classical sine-Gordon equation: solitons, which 
are time-independent and topologically non-trivial, and breathers, which are time-dependent and topo- 
logically trivial. Intuitively, a breather can be thought of as a bound state of a kink and an antikink 



,-m - Asin( -$ 
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oscillating in and out (i.e. breathing). Here, we will discuss only the solitons; the analysis for the 
breathers follows analogously. 

A major triumph of t he theory of non-linear partial differential equations was the construction of 
explicit solutions of (14E) for any number of moving solitons. The solitons' energies and momenta 
are conveniently expressed in terms of rapidities ctj , defined by Ej — M cosh aj and Pj — M sinh ccj , 
M = A 1 / 2 . The velocity of each is thus given by tanhaj (positive for a right-moving soliton). We have 
set the speed of "light" to be 1. 



Consider now a two-soliton solution of (149) on (—00, 00). This solution is usually expressed as 

'Xm(r) 



<fr(x,t) = 4 arg(r) = 4 arctan 



TZe(r) 



(150) 



where the r-function solution is given by: 

T = 1 



ei£2 (tanh^i^)" e -E 1 (x- a )-E 2 ( x -b)+P 1 t+P 2 t 

i \e ie -Ei{x-a)+P t t + e2e -E 2 {x-b)+P 2 t^ 



(151) 



The constants a and b represent the initial positions of the two solitons, and ej = +1 if the j soliton is 
a kink, while ej — — 1 if it is an anti-kink. 

What happens if we try to take the massless limit of this solution? For a wavepacket to have finite 
energy in the massless limit m — > 0, the rapidity \a\ must go to infinity. We thus define a = A + 9, 
and let A — > 00 such that the parameter m = \Me A remains finite. The energy and momentum of a 
right-moving "massless" soliton then reads 

E = P = me e (152) 
For a left mover, a = —A + 9, and its energy and momentum read 

E = -P = me~ e . (153) 
Suppose that both of these solitons are right-moving. Then the massless limit yields: 



l- eie2 e- A e- E ^- tt )- E ^- 6 ) +i 



-Ei(n-a) 



£2e 



-E 2 ( v -b) 



(154) 



where rj = (x — t) and 



A= -log^tanh^x-^)) 2 ] . 
This leads to an a priori strangely complicated solution of the Klein-Gordon equation. Observe that: 



arg{\- e l£2 e^v-^-E^-b) + ■ [ eie -^(n-a) + t ^-E^-h) j j 



arg 



1 + ie\e 



-Etiri-a) 



arg 



1 + if^e 



-E 2 (ri-b) 



— tan 



e x e 



-£1(77-0) 



tan 



£2e 



-E 2 ( v -b) 



(155) 



This is easily checked to be the sum of two one-soliton solutions; the factor A thus measures the extent 
to which the two-soliton solution is not a superposition of one-soliton solutions. 

More precisely, consider the limit a — > 00, 77 — > 00 so that Ei(rj — a) is finite. This corresponds to 
moving the first kink off to x = +00 and following it. The r function collapses to the one-kink form 
t = 1 + ic\e~ El ( :n ~ a ' . Moving this soliton through the second one corresponds to taking it to x = -co, 
or taking the limit a — > —oo,r/^ —00 (with £1(77 — a) finite). Discarding an overall multiplicative factor 
(which is irrelevant in the computation of <f> = 4arg(r)), we see that in this limit, r = l + ieie -Bl (' ,-a )~ A . 
Thus these preferred Klein-Gordon wave packets exhibit non-trivial monodromy. The foregoing time delay 
A is precisely the classical form of a massless scattering matrix. 

One obtains the same A for two left-moving solitons. For a left-moving and a right-moving soliton 
colliding one easily sees that the massless limit of (tanh 5(0(1 —012)) is unity. The solution collapses to 
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the superposition of a left-moving wave packet and a right- moving wave packet exactly as in ( 155 ), with 
no time delay. This is the classical manifestation of the fact that the left-right quantum scattering matrix 
Slr elements are at most rapidity-independent phase shifts. 

Consider now the Klein-Gordon equation on [— oo, 0] with the boundary condition ( 148 ). Here is 
a direct way of seeing the integrability of the Klein-Gordon (and indeed, sine-Gordon) equation with 
boundary conditions (148). The idea is to show that the method of i mage s can be used on (— oo, oo) even 
in the non-linear system, so as to replicate the boundary conditions (148) on [— oo,0]. The scattering of 
a kink, or anti-kink, from the boundary can be described by a three-soliton solution on (— oo, oo). These 
three solitons consist of the incoming soliton, its mirror image with equal but opposite velocity, and a 
stationary soliton at the origin, to adjust the boundary conditions (see (4^] for more details). If one takes 
the infinite rapidity limit of this three-soliton solution then the stationary soliton simply collapses to an 
overall shift of by a constant, while the mirror images (since they are moving in opposite directions) 
reduce to a superposition of two wave packets. One thus obtains: 



<f> = $o + 4 arg[\ + ieie 



-m- 



arg[l 



-E(r,-b)-\ 



where £ = x + t, rj = x — t. By direct computation one finds that this solution satifies (148) with 



-£l£2 



(2E + A) 
(2E - A) 



(156) 



(157) 



The constant Ab = —E(a + b) represents the delay of the reflected pulse. If one defines the classical 
boundary scale 9b via g = 2me~ SB , then this delay may be written as 



log 



-eje 2 tanh~(0 - 6 B ) 



(158) 



Note that the sign e 2 is to be chosen so as to make the argument of the logarithm real. This determines 
whether the reflection of a kink will be a kink or an anti-kink. Thus we see that 9b is the scale at 
which behavior crosses over from the region of the Neumann critical point (where the classical boundary 
scattering is completely off-diagonal) to the Dirichlct boundary critical point (where classical boundary 
scattering is diagonal). 

The conclusion of this section is clear: there are classical wave packets that scatter very nicely at the 
boundary. The price to pay is that they are considerably more complex than plane waves, and "weakly 
interacting" - that is, they scatter non trivially through one another. 

The natural attitude after having established such results classically is to see whether they are pre- 
served quantum mechanically. One can for instance establish integrability order by order in a loop 
expansion (which here amounts to an expansion in powers of (3). Here I want to show a more direct way 
to proceed, that generalizes to all sorts of theories. 



5 Quantum integrability 

5.1 Conformal perturbation theory 

As in the previous section, we start by considering theories with a bulk interaction. All what follows is 
based on the very insightful work of A.B. Zamolodchikov ptl|. Consider therefore the usual sine-Gordon 
model with actionp 1 ] 



l -Jdxdy {d x <S>f + (d y $) 2 + A J dxdy cos /?$ (159) 



*=2. 

As in the classical case, quantum integrability is established by proving the existence of non trivial inte- 
grals of motion. Rather than doing perturbation around (3 = 0, what one can do instead is perturbation 

21 It is safer to assume here that < |, so no counter term is necessary to define the perturbed action. Only a finite 
number of such counter terms would be required in general anyway, because of the "super-renormalizability" common to 
most perturbed conformal field theories. 
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around the conformal limit. This requires first the understanding that a conformal theory is integrablc 
P 2 ] , a rather straightforward property. 

Indeed, let us try to build a set of conserved quantities for a quantum field theory. We consider 
Euclidian space with imaginary time in the y - direction: a quantity will thus be conserved if its integral 
along two horizontal contours at different values of y gives the same result. Using complex coordinates, 
this will occur if we have a pair of quantities, say T n and 0„ such that dT n — <90„. Right at the conformal 
point, by analyticity, T, all its derivatives and (regularized) powers, do provide conserved quantities. 

To clarify this a little, let us consider the classical case. Going to imaginary time and complex 
coordinates, the equation of motion is 

d<9$ = ^ sin $ (160) 
Exercise: Show that the first pairs leading to conserved quantities are 

T 2 = (d$) 2 , 2 = -|cos$ 
T 4 = (9 2 $) 2 - i (c>$) 4 , 6 4 = j (3$) 2 cos $, (161) 

that is, the relation dTm — <90 2 n holds. Go back to real time, and find out which quantity is, indeed, 
conserved by time evolution. 

Away from the conformal point, what will happen is that there will sometimes be a deformation of 
(some of) these quantities that is still conserved. To see that, let us start by looking at the stress energy 
tensor, and see what dT becomes with a perturbation. To make sense of this question, we have to insert 
T inside a correlator, as usual. The difference with the conformal case is that now the action reads, quite 
generally 

S = S e ft + A J dxdyO (162) 

so we have to expand the Boltzmann weight in powers of A. This gives an infinity of terms, each of 
which has now a Boltzmann weight with a pure S c ft, so the results right at the conformal point can be 
used for them. But then it seems that T being analytic at the conformal point, nothing will make it non 
analytic away from it! That is not true because of what happens at coincident points. The integral of the 
perturbing field will affect only the behaviour near z, so we can use the OPE of T with the perturbation 

hO(z',z') dO(z>,z>) 
T{z)0iz , Z ) = j—^ + ——r- + ... (163) 

Now, using the identity ( fL6| ) it follows that 

dT = tt(1 - h)AdO (164) 

Hence, for T 2 = T, 9 2 = — 7rA(l — h)0, and we have a conserved quantity - to first order in A that is. 
Before wondering about higher orders, let us stress what ensured the existence of a conserved quantity: 
the fact that the residue of the simple pole of the OPE of T with the perturbation was a total derivative. 
One can then try to see whether there are other quantities for which a similar thing holds. Let us consider 
therefore the combination T4 = 4tt 2 : (<9</>) 2 : +A : (d 2 (f)) :. After laborious computation, one finds that 
the residue of the simple pole, in the sine-Gordon case of interest, is 

\ ZTT 577/ 

which is a total derivative when (notice this would still hold with (3 subsituted with ^) 

)-■>-( 3-^!-^ 



2 At least partly integrable - this subtlety does not seem to matter for most problem occuring in condensed matter. 
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The same sort of argument can be built to show that a conserved quantity can be obtained for every even 
n, which is the integral of a local field of dimension n. 

We now have to discuss what happens beyond first order. Suppose we carry out the computation to 
order n; a priori, we expect the result to be something like 

dT 4 = d<d 4 + 0(A n ) (165) 

where 84 is of order one in A. The left hand side has dimensions (4,1), while A n has dimensions 

(n(l — ^:), n(l — §^)j ■ This means that a local field with dimensions ^4 — n + n§^, I — n + n^J has 

to appear multiplying the A™ term. Since /3 2 < 87r, only a finite number of cases allow a positive set 
of dimensions, and for each of these except n — 1, and for 1 generic, one checks that there is no field 
with these dimensions. Hence, the conservation at lowest order extends, generically, to conservation at 
arbitrary order This proves quantum integrability, perturbatively that is. 

Finally, these conserved quantities also turn out to be in involution, ie they define mutually commuting 
operators. The proof is based on the Jacobi identity, which implies that, if two conserved quantities do 
not commute, then their commutator is also a conserved quantity - a little more thinking then establishes 
this is not possible. I also would like to remark that in the early literature about quantum integrable 
models, there is the implicit "suspicion" that quantization might destroy integrability, ie reduce the 
classical symmetry of the theory. It is important to realize that such a thing is not always true. In fact, 
in two dimensions at least, the quantum theory often has more symmetry than the classical theory. 

5.2 S-matrices 

The next step in the analysis of a massive integrable quantum field theory requires going to a scattering 
description |43[ . Let us assume quite generally that we have massive particles distinguished by some label 
a, with mass M a . We write their two momentum in terms of a rapidity variable a: E = M cosh a, P = 
M sinh a. 

We now consider scattering processes. There are the "in-states" , corresponding physically to a 
bunch of particles arranged on the x-axis by decreasing order of rapidities, which we describe for- 
mally by a ket \a±, . . . , ajv)S''"' a ™- At large time, they give rise to "out-states", formally described 

by \a'i, . . . , a4')o«t " ,a "', with an a priori different set of particles arranged on the x-axis by increasing 
order of, a priori different, rapidities (see figure 11). 



Pi p2 




Figure 11: Scattering of quasiparticles in an integrable, 1 + 1 quantum field theory. 

23 Non generic cases could be more complicated; consider the example of h = 1 for instance, where now arbitrary orders 
are allowed in the right hand side of (16q)! 
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Both the in and the out states are a complete set of states in a local quantum field theory, and they 
are connected by the S-matrix. 

Now, the existence of infinitely many conserved quantities has very drastic consequences on the 
scattering of these particles. Indeed, the conserved quantities have to act simply on the multiparticles 
states - they are in fact proportional to the sums of odd powers of the momenutm. As a result, it follows 
that, in the scattering process: 

(i) The number of particles is conserved; in fact, the number of particles of the same mass is conserved 

(ii) The final set of two-momenta coincides with the initial set of two-momenta 

From this in turn, it follows that the S matrix factorizes into a product of 2-body scattering processes. 
To see this, consider for instance conjugating the S matrix by the operator e ip i Xfi . Since the S matrix 
conserves momenta, it actually does commute with this operator, so we don't change anything. On the 
other hand, the operator has a non trivial physical action: it changes the space time coordinates of 
particle 1. If we chose p\ appropriately, we can arrange for particle 1 to scatter with the other particles 
only after they are very highly separated, so the scattering of this particle is a succession of two particle 
scatterings. By proceeding inductively, one deduces that indeed, the S matrix factorizes. Moreover, for 
the whole thing to be consistent, the scattering must be "associative"; that is, the scattering of three 
particles can be decomposed into three pairwise scatterings, with a result independent of which particular 
decomposition is used. This is illustrated graphically in figure 12 



Figure 12: Factorization of the scattering. 

This constraint, the so called Yang Baxter equation, is the pillar of the algebraic approach to integrable 
quantum field theories (for a review see for instance 

Notice that the YB equation is trivial if all the particles have different masses. It becomes more 
interesting in the case where several particles have the same mass, but differ by some other quantum 
number, eg the charge. In general, we define the S matrix elements by the relation (figure 13) 

I \in rya* ,a 2 I \out ?irr\ 

m,a 2 ) aua2 = SaU 2 2 |ai,a 2 ) / )0 , (166) 

Relativistic invariance contrains the S matrix to depend on the difference of rapidities ot\ — a 2 - 

In the case of the sine-Gordon model, perturbation in /3, considerations of (quantum affine) symmetry 

minimality and consistency assumptions, lead to the following results. The spectrum is made up of the 

kink and antikink of mass M oc A 2 - 2 s , together with breathers. For n — 1 < ~ < n, there are n — 2 such 

breather states. Their masses are Mk = 2Msin [kirg/2(l — g)]. 

The kink S-matrix is closely related to the matrix of Boltzmann weights in the 6- vertex model. There 

are three key amplitudes 

a(a) — sin[7(7r + ia)]Z(a) 
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a' 2 



a 1 
tX 1 



a' 1 




a 2 
<X 2 



Figure 13: The matrix element Sal'a 2 2 corresponds to the process illustrated here. 

b(a) — — sm(i"/a)Z(a) 
c(a) = sin(77r)Z(a), 



(167) 



where 7 = - — 1. 






Figure 14: The three types of possible processes involving kink and antikink. 



The element a (a 1 — 0:2) describes the process |«ia2) ++ —> \ a i a 2) ++ i as well as laio^) — * |aia:2) , 

b describes H > + — , c describes the non-diagonal process H ► — h (see figure 14), and there is a 

symmetry under interchange of kink to antikink (corresponding to $ — * — <f>). 



Exercise: Show that (167) gives rise, indeed, to a solution of the Yang Baxter equation. Hint: proceed 
graphically as sketched in figure 15 - to represent the matrix multiplications of YB, simply draw all the 
physical processes that connect a given pair of initial and final states, and add up their amplitudes. 

The function Z(a) is a normalization factor, which can be written as 



Z(a) 



1 



sin[7(7r + ia)] 



exp 



dy . taiy sinh[(7 - l)y] \ 
-00 2 V T sinhycosh(7y)y 



The breather-kink and breather-breather S matrices are well known ; we do not write them down here 
(see below for some examples, and in the appendix). 

These S matrices of course are not only characterized by the fact that they must solve the Yang- 
Baxter equation: there are several other physical requirements, like unitarity and crossing symmetry. In 
the formula above, these translate into the relations 



a(a) = b(in — a) 
c(a) — c(iTr — a) 



(168) 



form which unitarity S(a)S(—a) follows. 

In addition, one must have a "closed bootstrap" : for instance, breathers appear as poles in the kink 
antikink scattering, and their S matrix can be computed using this fact (for a review see @). 

When 1/g is an integer, the bulk scattering is diagonal (c vanishes) and a — ±6. Therefore, the 
only allowed processes are transmissions: particles go through one another without exchanging quantum 
numbers. This is the simplest case, to which we will restrict in what follows. 
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<X 3 



<X 2 




Figure 15: The three types of possible processes involving kink and antikink. 



It is reasonable to think of the S matrix, together with the mass spectrum, as the quantum equivalent 
of the knowledge of the r function in the classical case 0. We can then take the massless limit exactly 
as we did before, by letting the mass parameter M — > 0, and at the same time boosting the rapidities. 
One obtains then a collection of left and right moving kinks, antikinks and breathers. The LL and RR 
scattering matrices take exactly the same expressions as in the massive case in terms of the new rapidities 
9, while the LR S matrices go to constants, which can in most cases just be forgotten. This provides an 
alternative description of the free boson in terms of "massless scattering" . We will see in a little while 
how quantities of the conformal field theory can be recovered, if one wishes, within that description. 

It is fair to stress here that the RR or LL scattering are hard to make sense of in the context of a true, 
physical scattering process. If particles have the label R say, this means they are moving at the speed 
of light in the right direction, so, for instance, all the lines in the figure 11 which illustrates the Yang 
Baxter equation become parallel! The point is that the scattering has to be interpreted in the massless 
limit as a set of commutation relations for creation operators, an idea which we will discuss below. The 
idea of massless scattering appeared a bit weird at the beginning jl7| , (|o) , [Q , but its predictive power 
and favorable comparison with experiments gained it respectability quickly. A more rigorous approach 
along the lines of lattice regularizations is proposed in [49) . 



5.3 Back to the boundary sine-Gordon model 

For the moment, we finish following the logic of the classical analysis in the quantum case. First, one 
can prove that the quantum sine-Gordon model with a bulk and boundary interaction 



S = 



1 



<// /./,-/ (<9 X $) 2 + {dy^f + Acos/3$(x,?/) +A / dycos^$(0,y). (169) 



allows the existence of conserved quantities as well |H] . The proof proceeds in the same spirit as for the 
bulk case, and we are not going to reproduce it here, though it is an excellent exercise for the diligent 
reader. 



24 There are more accurate formulations of this statement. 
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Exercise: Show that, in the classical theory, the pair T4, 02 still leads to a conserved quantity provided 
there exists still another local quantiy #3 such that 



r4 + e 2 -r 4 -e a | ai=0 = -< 



(170) 



Show then, that for the sine-Gordon model with a potential = 0, y)] at the boundary, only the 

choice V cx cos + $0) allows the existence of the integral of motion. 

Integrability means that the particles also have to scatter nicely at the boundary, ie scatter one by 
one, without particle production, in a way that is compatible with the bulk scattering. This latter 
condition can simply be expressed graphically as shown on figure 16, and corresponds to the "boundary 
Yang-Baxter" equation. 





Figure 16: The boundary Yang-Baxter equation. 



The problem of determining the reflection matrix is purely technical, 
answer in the massless limit. Introduce 



It does have quite a simple 



R+-(6) =R-+(0) 



-i exp(76>) 



exp[ixg(6)] 



1 — i exp(7#) 

R ++ (6) =R—(0) = . 1 , m expfa g (0)] 

1 — iexp(7#) 

where \g is a phase that will disappear at the end of the computations, 7 = - — 1. The reflection matrix 
for coupling A is then given by R(9 — 9 b) where Tb = e 6s cx A 1 '^ 1- ^ (the exact correspondence depends 
on the regularization scheme; it is given in |]36f in the case of dimensional regularization) is the equivalent 
of the Kondo temperature. As 9 — > 00, R++ — > (the scattering is completely off-diagonal) corresponding 

to Neumann boundary conditions, while as 9 — > —00, R^ — > (the scattering is completely diagonal) 

corresponding to Dirichlet boundary conditions. Notice also the unitarity condition \R^ — | 2 + \R++ 1 2 = 1. 



6 The thermodynamic Bethe-ansatz: the gas of particles with 
"Yang-Baxter statistics" . 

The thermodynamic Bethe ansatz was probably written first in [j52| in the context of the XXZ model. 
Its use in quantum field theory, in particular to compute the central charge and study RG flows, was 
pioneered by a beautiful series of papers of Al. Zamolodchikov; see for instance p3|]. A useful and 
pedagogical review on many of these topics can be found in 
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6.1 Zamolodchikov Fateev algebra 

It is convenient to think of the particles in terms of creation and annihilation operators. For this, let us 
introduce, still denoting the type of the particles by a label a, operators Z a (0) and Z\(0) satisfying the 
relations 

Z ai (9l)Z a2 (9 2 ) — S ai a 2 (01 — 02)Za 2 (#2) Z ai (01 ) 

< (01)4,(02) = S aia2 (0i -9 2 )Zl 2 (0 2 )Zl(9i) 

Za^zl^) = s aia2 (e 1 -e 2 )zl 2 (e 2 )z ai (e 1 ) + 27r6 aia2 6(e 1 -e 2 ) (171) 

Here, we restricted to the case of diagonal scattering. Note that the compatibility between the first two 
relations uses unitarity in the form S'(9) = S~ 1 (9) = S(—9). The space of states is generated by the kets 

\9 1 ,...,0 n ) au ^ an =Zl(9 1 )...Zl(9 n )\O), (172) 

where |0) denotes the physical vacuum. Similarly, the dual space is generated by the bras 

(0„,...,0i| = (0\Z a J8 n )...Z ai {8 1 ) (173) 



The metric is, from (171), induced by 

Ol (0i|02> O2 =2mJ aia2 <J(0i- 2 ) (174) 
If for instance 9\ > 9 2l then the in and out states are, respectively 



|0i,02)a" a2 - \9\,9 2 ) ai 



a 2 



|0i ! 02CL = |02 ! 0i) OlO2 (175) 

When the rapidity sets are not ordered, one obtains states which are neither in nor out; or course they 
are related to either of these by products of S matrix elements. 

To make things more concrete, let us discuss briefly wave functions in coordinate representation, 



restricting for simplicity to two particles. To satisfy the relations (171), it is easy to see that the wave 



function must have a singularity at coincident coordinates, and be of the form 

\W ai a 2 <* I dx X dx 2 e «**i+ft*») 13^)^-^(0! -0 2 ) / d Xl dx 2 e «**i+^> \ Xl ,X 2 ) aia2 

J x\<x 2 J X\>X 2 

(176) 

where we assumed that the particles are fermions, S(0) = — 1. Equivalently, one has 

' e ^x l+ p 2 x 2) + SaiaAQi _ ^+^1] |0i,0 2 )^ O2 (177) 

where we see the appearance of the well known Bethe wave function pa] . 



d9 x d9 2 

>e 2 



6.2 The TBA 

The next step is to get a handle on the massless scattering description. The latter turns out to be quite 
convenient to discuss thermodynamic properties, and this is what we shall start with. 

As a simple example we consider a hypothetical theory made up of a single type of massless particle, 



say right-moving, with energy and momentum parametrized as in (152). The scattering is described by 
a single S'-matrix element Sjm. Quantizing a gas of such particles on a circle of length L requires the 
momentum of the ith particle to obey (we have set h = 1) 

exp (ime^L) ]J S RR (9 t - 9j) = 1. (178) 
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One can think of this intuitively as bringing the particle around the world through the other particles; 
one obtains a product of two-particle S-matrix elements because the scattering is factorizable. A bit 
more rigorously, one can deduce this from the wave function as in ( 176 ) 

Going to the L — > oo limit, we introduce the density of rapidities indeed occupied by part icles p{9) 
and the density of holes p. A hole is a state which is allowed by the quantization condition ( |178] ) but 
which is n ot oc cupied, so that the density of possible rapidities is p{6) + p h {9). Taking the derivative of 
the log of (Il7|) yields 



2tt[p(0) + p h {0)] = mLe 



K(6-0')p(6')d6', 



(179) 



where 



To determine which fraction of the levels is occupied we do the thermodynamics, following the pioneering 
work of Yang and Yang. The energy is 



and the entropy is 



S 



p(6)me e de, 



[(p + p h ) \n(p + p h ) - pHp) - p h ln(/)] M. 



Exercise: derive this relation by using Stirling's formula T(z) 



2 e 



2tt. 



The free energy J- = {£ — TS) is found by minimizing it with respect to p. The variations of £ and 
S are 

5£ = X!^ 8pme e d9 

SS = f™^ [(5p + 5p h ) ln{p + p h ) - Sp ln(p) - Sp h ln{p h )] dO 

It is convenient to parametrize 
giving 



(180) 



<5pln(l + e £ / T ) +8p h Ya{l + e-*' T 



dO. 



Using (|179J) allows us to find p in terms of p. Denoting convolution by *, this gives 2it(8p + 8p ) = K-kSp 
so 



SS 



Hence the extremum of T occurs for 



e K 
T + 2^ 



★ ln( 1 + e~ e/T 



Spde. 



me a = e + T-- * In ( 1 + e 



e/T 



and one has then, expressing p h from (|l79| ) and using (181) 



T 



-LT 



2ttT 



e" In 1 



£ / T ) dO. 



(181) 



(182) 



It is a simple exercise to show that this formula, together with (181), generalizes to a theory with 

several species of particles, provided the scattering is diagonal. This corresponds to the case = g = j, 
t an integer, to which we restrict in what follows. In that case, recall that we have a kink and antikink of 
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mass parameter m, and breathers of mass parameter = 2m sin 2 (t-i) > w ith k = 1, • • ■ ,t — 2. We will 
also allow for different chemical potentials fik for the various particles. Defining now e's through 

" W -«pfc^l (183) 



the equivalent of ( |179| ) is now 

/oo 
K jk (9-6')p k (6 , )d9', (184) 
-oo 



and the equivalent of (181) 



rune = e 



The equivalent of ( 182 ) is, in turn 



g = £-TS-Y J ^M k = -LT 2 Y J Pk r e e ln(l + e^)d9. (186) 
k k Znl J -°° J 



For the case g = \ for instance, one has 

K bb = 2K ++ = 2K, 



cosh( 



z— cosh 8 

^ + =K +h = -2^^- Q (187) 

It has become common in the literature to reformulate the TBA in a convenient form by using simple 
diagrams. It is a laborious but straightforward exercise to demonstrate, using the kernels given in the 
appendix, that (185) is equivalent to the following simple system ^| 



* = TJ2 N ^Y * ln (l + e^) (188) 



Here, s(6) = CO gh(t-i) W' = 1 ^ t ne nodes j and fe are neighbours on the following diagram, 
otherwise 

P + 

1 2 s t-3 / 

O 0-- --0---0 O t-2 



Exercise: establish this for the case g = \- 

The equations ( |188| ) have to be supplemented by the boundary conditions 

ej re m je e , 6» » 1 (189) 



25 The case where g is not of the simple form 1/ 'integer can also be handled of course. It is technically more difficult 
because the scattering is non diagonal, so an additional Bethe ansatz is necessary to diagonalize the scattering to start with, 



before the periodicity of the wave function can be imposed j54j . 
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6.3 A standard computation: the central charge 

The thermodynamics of a chiral theory like the one we just studied is not so exciting; this is because, after all, the 1 + 1 
theory is conformal invariant, so the results at different temperatures are essentially equivalent. This can easily be seen on 
the TBA equations: a change in T, or the mass scale m, can be fully absorbed by a boost of the particles, ie a shi ft of 
rapidities, exactly like for the changes in mass scale encountered before. As a result, we see that the integrals in (182j) are 
independent of the temperature, so T oc LT 2 . 

That this is so, and the coefficent of proportionality, are directly related with considerations from the beginning of these 
lectures. Indeed, we have T = — TlnZ, where Z is the partition function of the one dimensional quantum field theory at 
temperature T. In Euclidean formalism, this corresponds to a theory on a torus with finite size in time direction R = X/T. 
By modular invariance, identical results should be obtained if one quantizes the theory with R as the space coordinate 
. For large L, Z = e ~ E ( R ) L ^ where E(R) is the ground-state (Casimir) energy with space a circle of length R. Thus 
T = LE(R) I R. Conformal invariance requires that at a fixed point this Casimir energy is E(R) 

central charge. 



Going back to the thermal point of view, T = — Hf ^jr- and the specific heat isp^ C ■ 



-■5s, where c is the 

6-R ' 



It is possible to analytically find this central charge from (181). T his i s a bit technical, but worth studying 
crucial a, po steriori test of the whole thing. We take the derivative of (181) with respect to 9 and solve for e? . 
this in (182), we have 



, since it is a 
Substituting 



^ = I <» [§ Wl + e-«/r) - / g ln(l + e-W/*W ~ & TT^JJT 

= - 5£ / [mi + e - e/T ) + ^ i^^W 



1 dx_ 

dO 



(190) 



where we use (181) again to get to the second line. We can replace the integral over 9 with one over e, giving an ordinary 
integral 



TL 
Air 



ck 



(-oc) 



ln(l + e~^ T ) + 



e/T 



l + e e W/ T 



A change of variables gives 



1 



T 2 L 

T= C , 

2n \ 1 + xo 



(191) 



where C(x) is the Rogers dilogarithm function 

^ ) = -if( 



+ I dy, 



1-2/ 



and xq = exp[e(— oo)/T] is obtained from (181) as 



- = ( 1 + -)'" 

XQ \ XfjJ 



with I = ± [ K. 

For example, when the S matrix is a constant, K = 0, xq = 1 and 



T - 



LT 2 



21 



(192) 



(193) 



where we used L(l/2) = 2— Here we find Ci = \. In a left-right-symmetric quantum field theory, the right sector makes 
the same contribution, giving the total central charge c=\ required for free fermions. 

Similar computations can be carried out for more complicated theories, leading to beautiful expressions of central 
charges in terms of sums of dilogarithms (see eg |56|). In the case of interest, one finds of course c = 1. 



6.4 Thermodynamics of the flow between N and D fixed points 

We now wish to do the thermodynamics in the presence of the boundary, to obtain the boundary free 
energy, and the associated flow of boundary entropies. To start, it is better to map the problem onto 
a line of length 2L {—L < x < L) by considering the left movers to be right movers with x > 0. Thus 
we have only R movers scattering among themselves and off the boundary, which can now be thought 
of as an impurity (a particle with rapidity Ob)- The reflection matrix becomes a transmission matrix, 
with appropriate relabellings, for instance R-\ — > etc. (This trick is the same than what we did 

26 With massive particles or with nontrivial left-right massless scattering, T does depend on MIT, giving a running central 
charge. 
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for boundary conformal field theory, and can only be used in the massless limit.). For simplicity, we put 
periodic boundary conditions on the system; these do not change the boundary effects at x = 0. 

Recall we consider only the case 7 = - — 1 a positive integer, where the bulk scattering is diagonal. 



The impurity scattering still is not, but we can redefine our states to be |1, 2) 
the impurity scattering is now diagonal: 



Tim 

T 2 2{0) 



R++ 
R ++ 



R+- 
R+- 



exp [ixg(0)] 



tanh 



2 



exp [ixg(0)} 



± \-))/V2 so that 



(194) 



We can now write the Bethe equations. These differ from the bulk ones only by the presence of the 
additional impurity scattering 



2n( Pj (9) + ${B)) = mj e e + £ K jk *p k (6) + ± Kj {6 - 9 B ) 



(195) 



where 



K 



Kj(6-d B ) 



1 d 

~i~dB 
1 d 

id§ 



(196) 



The effect of the boundary is seen in the last piece of (195) proportional to 1/L 



The minimization equations are independent of the boundary terms, since these do not appear directly 



in £ or S, and they disappear when one takes a variation of (195). Thus equations fll85j) still hold 



Boundary terms do enter the free energy or the grand potential however when one rewrites it in terms 
of the e's. One finds 



T = T } 



bulk 



T 



de 



2tt ^ 

3=1 



ln(T/T B ))ln(l + e- 6 'W). 



(197) 



As discussed before, J-'bulk = — ^§-T 2 L in a massless bulk theory, where c is the central charge of the 
conformal field theory, c = 1 here. The second term in (197) is the boundary free energy. 

Although the equations ( |185| ) for e(9) cannot be solved explicitly for all temperatures, the free energy 
is easy to evaluate as T — > and T — > 00, as we will show next. Moreover, one can extract the analytic 
values of critical exponents by looking at the form of the expansions around these fixed points. Also, 
they are straightforward to solve numerically for any T. 

Several notes of caution are necessary. At the order we are working, the formula for the entropy is not 
quite correct, because there are 1/L corrections to the Stirling formula used in its derivation. Also, at 
this order, the logarithm of the partition function is not £ — TS: it depends not only on the saddle point 
value of the sum over all states, but also on fluctuations. Their net effect is that we cannot compute the 
g factors from T alone. However, both of these corrections are subleading contributions to the bulk free 
energy, and do not depend on the boundary conditions. Therefore we can still compute differences of g 
factors from T\ the corrections are independent of the boundary scale 9b and cancel out of the difference. 

We can evaluate the impurity free energy explicitly in seve ral limits. In the IR limit T/T B — * the 
integral is dominated by 9 — > 00 where the source terms in ( 185 ) become very big. Hence e r (oo) = 00 and 
the impurity free energy vanishes in this limit. In the UV limit T/Tb — > 00 the integrals are dominated 
by the region where —9 is large so that the so urce terms disappear in (185) and the e r go to constants. 
These are found by using the alternative form (|188|), which reads here, denoting Xj = e Cj / T 



(198) 



k 



27 If 7 is even, this actually makes the bulk scattering completely off-diagonal (e.g. |11) scatters to |22)), but the TBA 
equations turn out the same. 
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One finds 

x n = e e ^-°°V T = (n + 1) 2 - 1; 11,2 = 7 (199) 

Therefore we obtain 

-Fir, 



9D 



uv 



IR 



= E^=i-f (n) M 1 + l/x„) + )ln(l + l/a; ± ) (200) 

where 

J(r) = / g Kr(0 ) = *((,). 

From results given in the appendix, one finds JO) = n/2 and /(+) + = 7/2, and thus 

9n _ 7, 7+ 1 ,V^ n i ("+ 1 ) 2 



n„ 9 m A^l 9 



on 2 7 ^2 nfn + 2) 

n— 1 7 

= iln( 7 +l) = ilnt. (201) 

This is in agreement with the ratio calculated from conformal field theory. 

We can also find the dimension of the perturbing operators. From the equations one deduces the 
following expansions for T/Tb large: 

Y r (ff) = e £ " (9) = ^y r li) e -W/(7+i], 
3 

As a result it is straightforward to see that near A = 0, T can be expanded in powers of (Tb/T) 27 ' 
On the other hand we expect T to be an analytic function of A 2 . Hence 

Aoc(me e -) 7/(7+1) . (202) 

This agrees with the conformal result that the perturbing operator cos[/3$(0)/2] has boundary dimension 
d = 1/(7 + 1) = P 2 /8tt. In the IR limit of T/Tb small, one can expand out the kernels k t in powers of 
exp(6*b — 6). This leads to the fact that the irrelevant operator which perturbs the Dirichlet boundary 
conditions has dimension d=2. This is the energy-momentum tensor. (Recall that there is another 
irrelevant operator in the spectrum with dimension d = 7 + 1, which for < 7 < 1 is the appropriate 
perturbing operator.) 



7 Using the TBA to compute static transport properties 

Let us pause for a moment to compare the gas of Yang-Baxter interacting quasi particles to say free 
fermions. Within the TBA, the interactions have been fully encoded into non trivial pseudo energies 
e.j{0): that is, at temperature T, the filling fractions of the various species are not independent, but 
correlated via the coupled integral equations discussed previously. This has some striking consequences. 



For instance, we see from (199) that the filling fraction of kinks or antikinks at rapidity —00 (ie at 
vanishing bare energy) is / = j. Except for t = 2 (which is a free fermion theory) there is no symmetry 
between particles and holes. It is important to realize that the interactions would have other effects, in 
general, for other questions asked. For instance, in the case of free fermions, the total density n = p + p h , 
p = nf, the fluctuations also depend on the ej through the well known formula (A/?) 2 = nf(l — /). Such 
a formula does not hold in the present case: the fluctuations of the various species are correlated - their 
computation plays an important role in the DC noise at non vanishing temperature and voltage, see (57). 
Similarly, physical operators have complicated matrix elements in the multiparticle basis; the current for 
instance is able to create any neutral configuration of quasiparticles by acting on the vacuum. There is 
thus a somewhat deceptive simplicity in what we have done so far. However, for the DC conductance, it 
turns out that the knowledge of the distribution functions is all that is necessary, so for that particular 
aspect, our quasiparticles are not so far from free ones. 
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7.1 Tunneling in the FQHE 

At this stage, it is useful to recall the tunneling problem of the introduction: we had L and R moving 
electrons that were backscattered by the impurity. Certainly if a R moving particle bounces back on 
the gate voltage to become a left one, the charge Qr + Ql is conserved. Now Ql + Qr. is essentially 
the charge of the even field in the manipulations discussed in the introduction, which we found has no 
dynamic indeed. Now when a R mover bounces back as a L mover, there is a change in the non conserved 
charge Q = Qr — Ql] this one is proportional to the charge of the odd field, which has a non trivial 
dynamics. Now, following carefully the formulas for bosonization, one finds the simple result that a right 
moving kink, for which ^ J d x $ = 1, also has physical charge Q = 1, and similarly for antikinks and 
left moving particles. Therefore, the non conservation of the physical charge due to backscattering is 
the same as the non conservation of charge in the boundary sine-Gordon model. More precisely, when 
a kink comes in and bounces back as an antikink, as happens most of the time near the UV fixed point 
(Neumann boundary conditions), the charge Q is conserved in the original problem. On the other hand, 
when a kink bounces back as a kink, as happens near the IR fixed point (Dirichlct boundary conditions), 
the charge in the original problem is not conserved; rather, AQ = —2. Let me stress here that the kink 
in the boundary sine-Gordon theory however would look horribly complicated in the original problem, 
because the changes of variables we have performed are non local. Only the conserved charge is easy to 
follow 



7.2 Conductance without impurity 



In the absence of impurity, that is with Neumann boundary conditions in the original boundary problem, 
charge is straightforwardly transported. A right moving kink or antikink just goes through. Of course, 
if there are as many particles of each specie, no current is transported overall. If however, a voltage V 
is applied, kink and antikink are at a different chemical potential, fi = ±-|- - this follows since the U(l) 
charge in the boundary sine-Gordon model is nothing but the physical charge Q. The current that flows 
through the system is thus 



/ = 



(p+-p-)(0)d0 



We can use our TBA to evaluate this expression quickly. First, we introduce the filling fractions 

1 



f± 



l + e (e±T%)/T 

Second, we observe that the very convenient identity rij = pj 



(203) 



(204) 



Pi 



1 dej 
2ir dB 



e+ = e_ = e. 

Exercise: Prove these two statements by staring at the TBA equations. 
It thus follows that 



1 



2tt 



dc 



and thus 



T 
2^ 



d0—\n 
d9 



dB 

e -V/2T e -e/T 



holds, and that, moreover, 



(205) 



„V/2T e -e/T 



(206) 



The current is thus entirely determined by the values of e at ±oo, exactly like for the central charge. As 
before, e(oo ) = oo, but the value of e(— oo) now does depend on the voltage. One finds in fact, solving 
again (188) but with a voltage, 



sinh(n + l)V/2tT 



sinh V/2tT 



1. 



e±(-oo)/T _ 



sinh(< - l)V/2tT 
sinh V/2tT 



(207) 



28 Charge is like current here, where we have set the Fermi velocity equal to one. 
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(observe one recovers the result ( 199 ) as V — > 0) from which an elementary computation shows the simple 
result (recall g = \) 

I= 9 — (208) 

The bizarre factor of 2n occurs here because we have set e = % = 1 (recall that in physical units, 
I = 9^V). 

Exercise: Prove the last two formulas. 

This is what one expected, and of course there are quicker ways to derive this result. The point 
however, is that the same computation carries over without much additional difficulty to the case where 
the impurity is present. 

7.3 Conductance with impurity 

In the general case, we will write the source drain current as I = Iq + Ib where Iq — is the current 
in the absence of backscattering, and Ib is the backscattered current. In the original problem, Ib is for 
instance the rate at which the charge of the right moving edge is depleted. Of course, d t Qh = —dtQn 
in each hopping event, so Ib = dt%- I n the steady state this rate is constant. When for instance V is 
positive, there are more kinks than antikinks injected with a thermal distribution into the system from 
their respective infinite reservoirs; it is assumed that these reservoirs are so big that the backscattering 
does not change their properties. 

We now derive an analytic expression for this backscattering current using a kinetic rate equation for 
quasiparticles of the Bethe ansatz. It is possible to compute the rate of change of AQ/2 in the basis of 
the Bethe ansatz quasiparticles, since each scattering event of a kink (antikink) into an antikink (kink) 
changes the physical charge AQ/2 by —1 (by (+1)). This kinetic equation is of course very familiar. 
However, in general there would be no reason why it should be applicable to an interacting system. But 
it is exact in the case that we are considering, even though the system is interacting. The reason for this 
lies in the constraints of integrability: as discussed above, in the very special quasiparticle basis of the 
Bethe ansatz, these quasiparticles scatter off of the point contact independently ("one-by-one"), and all 
quasiparticle production processes are absent p| . 

This allows us to express the rate of change in AQ, in terms of the transition probability |T_| | 2 

(recall that in the unfolded point of view = R++) and the number of kinks and antikinks (carriers 

of charge AQ = ±1) in the rapidity range between 9 and 9 + d9 

n±(9)f ± (9)d9, 

where n± is the density of states and f± are the filling fractions. The number of kinks of rapidity 9 that 
scatter into antikinks per unit time is 

|T + _(0)lV-^ (209) 

where p^ is the probability that the initial kink state is filled and the final antikink state is empty (in 

all these quantities, there is also a V dependence, which we keep implicit here). For a system of free 
fermions we would have 

P+ -=f+(l-f-) 

but in our interacting system we only have 

P+- = f+ - /+- 

where /_| is the probability that both, the kink and the antikink states are filled. For the number of 



antikinks of rapidity 9 that scatter into kinks per unit time one finds a formula similar to (209), with 
p_l — > p |_. In the final rate equation, only the difference between these two probabilities 

P+- - P-+ = f+ - /- 



29 Some more detailed justifications are available; see |58| and references therein. 
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appears. Notice that the unknown /_| = / |_ has cancelled out (it can in fact be determined by 

techniques more elaborate than the TBA [p7|). Therefore, the backscattering current is 



Ib(V) 



d9n(9)\T + ^(9-d B )\ 2 [f + (d)-f-( 



(210) 



All ingredients in this formula are exactly known: the scattering matrix has a simple analytic form and 
the occupation factors and densities of state are obtained exactly from the thermodynamic Bethe ansatz 
(TBA). Notice that this equation is valid for any value of the driving voltage V. It thus automatically 
describes non-equilibrium transport. 

By the same manipulations as before, it then follows that 



(16- 



i 



2(t-l)(0-e B ) d g 



In 



1 + e- y / 2T e- £ / T 



1 + e V/2T e -e/T 



(211) 



Of special interest is the linear conductance, which we obtain by taking a derivative at V — 0; this gives, 
after reinserting the factor 2-7T, 



G 



(t-1) 



d6 



1 



1 + e £ / T cosh 2 [(i - 1)((9 - ln(T B /T))] 



(212) 



The resulting curve is shown in figure 17, together with experimental results |13| and the results of Monte 
Carlo simulations |L4| , for g = | . 
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Figure 17: Comparison of the field theoretic result with MC simulations and experimental data for g 



The agreement with the simulations is clearly very good (the is one and only one fitting parameter - 
the horizontal scale - , accounting for the unknown, non universal ratio of the experimental gate voltage 
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(a "bare" quantity) to the parameter A in our renormalized field theory) . As far as the experimental data 
go, it is also very satisfactory, except in the strong backscattering regime. Recall however that the field 
theoretic prediction holds true only in the scaling limit: the experimental data are still quite scattered 
for low values of G, indicating that this limit is not reached yet - actually the "noise" is of the same order 
of magnitude as the discrepancy from the theoretical curve, as reasonably expected. 

Exercise: The problem had been solved previously in the simplest case of g = 5, where one can 
refermionize the hamiltonian for the boson (f> e . Look at this solution, and compare with what we have just 
done: what is the meaning of kink and antikink, what is the bulk scattering, the boundary scattering? 
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Conclusions: further reading and open prob- 
lems 



I am now leaving you at the beginning of a very exciting domain. Let me suggest some further reading 
and open problems. 

These lectures have stopped short of really tackling the problem of boundary fixed points classification. 
Equipped with what you learned here, you should not have much difficulty reading the paper of J. Cardy 
on boundary states This paves the way to questions that are still open. For instance, the problem 
we have studied at length generalizes, for tunneling in quantum wires where the spin of the electrons has 
to be taken into account, to a "double sine-Gordon problem", involving two bosons. Surprisingly it has 
been shown ||Io|| that new non trivial fixed points do exist in that case, besides the obvious Neumann and 
Dirichlet possibilities. With a few exceptions Q, nobody knows how these fixed points precisely 
look like! As I will mention again below, what we have discussed is also very close to the Kondo problem. 
You can learn more about fixed points and conformal invariance by reading the papers of Affleck and 
Ludwig on the multichannel Kondo problem fl6| . There, you will also discover an aspect that I have 
neglected by lack of space: how multipoint correlators can be evaluated at fixed points by further using 
conformal invariance [l8| . 

The integrable approach can also be pushed further to allow the computation of AC properties, 
together with space and time dependent Green functions, in the cross-over regime. The idea here still 
relies on massless scattering; but now, one has to evaluate matrix elements of physical operators, and 
these are usually pretty complicated. Moreover, an infinity of these matrix elements are a priori needed: 
for instance, the current operator is able to create any neutral configurations of quasiparticles out of the 
vacuum! It turns out however that, first, the matrix elements can be determined by algebraic techniques 
J59| , ]6C| ] (the latter reference is recomended as a first reading; the first is a bit hard to read), and second, 
in many cases, only a few of these matrix elements are required to obtain controlled accuracy all the 
way from the UV to the IR fixed point. Using that technique, for instance the current current correlator 
itself can be evaluated, at least at T = V = |3l| (there does remain a non trivial dependence upon 
space, time, and the coupling A). Determining correlators with a finite temperature or voltage is still 
more difficult; some progress in that direction has been made [j62j] . pSgjj , but a lot remains to be done. 

In another direction, for those of you who are more formally oriented, it should be clear that what I 
just described is the tip of an iceberg of beautiful mathematical structures: see ||^],(6^] and the series 
. Let me just mention here that the Kondo problem, which would be described by 

H = H e = ^J° dx [n 2 + (d x <S>) 2 ] + A (Ve-^°>/ 2 + S- e ^°)/ 2 ) , (213) 

is just around the corner: it actually does have deep relations with the boundary sine-Gordon model, and 
with the subject of quantum monodromy operators. Especially exciting results have actually appeared 
recently, concerning an exact duality between the UV and IR regimes of the problem l]68|| , f69| , and 
exhibiting tantalizing relations with the recent breakthroughs in 4D SUSY gauge theories |7C| . 

It is also fair to stress that the methods developed within the context of quantum impurity problems 
can be generalized to different systems of physical interest in 1 + 1 dimensions: an example is the amazing 
recent mapping of the two-ladder problem onto an 5*0(8) Gross Neveu model |n]]. It is very likely indeed 
that more such problems are awaiting us in the near future. 

Finally, the traditional question is, can any of this be generalized to more than 1 + 1? Well, the 
recent excitments in string theory are centered around somewhat similar ideas in 3 + 1, where, roughly, 
integrability is replaced by supersymmetry, an incredibly powerful tool. As for 2 + 1,1 don't quite think 
it's over yet. 
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Appendix: Kernels 

We use the following convention for Fourier transform: 

/°° Aft 

where 7 = t — 1. The bulk kernels Kjk arc well known; they can be written in the form (± stand for kink 
and antikink) 

K jk = S 2 cosh y cosh( 7 -j)ysmhfc g fc = 1 _ fc 

cosh 7?/ sinh y 

~ cosh y sinh jy 

K j,± = 



cosh 7j/ sinh y 
~ - sinh( 7 - l)y 

2 cosh jy smhy 

(215) 



with Kjk = Kkj. The boundary kernels are 

sinh j y 



2 sinh y cosh jy 



sinh(7 - l)y 1 



2 sinh 2y cosh jy 2 cosh y 

* + = o^l 7 " 1 ^ ■ ( 216 ) 
2 smh 2y cosh 72/ 

(217) 



Finally, s 



2 cosh y ' 
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